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ABSTRACT 

CAIRNS (Cluster And Infall Region Nearby Survey) is a spectroscopic survey of the infall regions 
surrounding nine nearby rich clusters of galaxies. In Paper I, we used redshifts within ^ 10/i~^Mpc of the 
centers of the clusters to determine the mass profiles of the clusters based on the phase space distribution 
of the galaxies. Here, we use 2MASS photometry and an additional 515 redshifts to investigate the 
environmental dependence of near-infrared mass-to-light ratios. In the virial regions, the halo occupation 
function is non-linear; the number of bright galaxies per halo increases more slowly than the mass of 
the halo. On larger scales, the light contained in galaxies is less clustered than the mass in rich clusters. 
Specifically, the mass-to-light ratio inside the virial radius is a factor of 1.8 ±0.3 larger than that outside 
the virial radius. This difference could result from changing fractions of baryonic to total matter or 
from variations in the efficiency of galaxy formation or disruption with environment. The average mass- 
to-light ratio M/Lk = 53 ± 5ft. implies 51„i = 0.18 ± 0.03 (statistical) using the luminosity density 
based on 2dFGRS data. These results are difficult to reconcile with independent methods which suggest 
higher Reconciling these values by invoking bias requires that the typical value of AI/Lk changes 
significantly at densities of ^Spc- 

Subject headings: dark matter — galaxies: clusters: individual (A119, A168, A194, A496, A539, A576, 
A1367, A1656(Coma), A2197, A2199) — galaxies: kinematics and dynamics — 
cosmology: observations — infrared:galaxies 



1. INTRODUCTION 

The relative distribution of matter and light in the uni- 
verse is one of the outstanding problems in astrophysics. 
Clusters of galaxies, the largest gravitationally relaxed ob- 
jects in the universe, are important probes of the distri- 
bution of mass and light. Zwicky (1933) first computed 
the mass-to-light ratio of the Coma cluster using the virial 
theorem and found that dark matter dominates the clus- 
ter mass. Recent determinations using the virial theo- 
rem yield mass-to- light ratios of M/Lsj ^ 25O/iMq/L0 
(Girardi et al. 2000, and references therein). Equating 
the mass-to-light ratio in clusters to the global value pro- 
vides an estimate of the mass density of the universe (Oort 
1958); this estimate is subject to significant systematic 
error introduced by differences in galaxy populations be- 
tween cluster cores and lower density regions (Carlberg 
et al. 1997; Girardi et al. 2000). Indeed, some numerical 
simulations suggest that cluster mass-to-light ratios ex- 
ceed the universal value (Diaferio 1999; Kravtsov & Klypin 
1999; Bahcall et al. 2000; Benson et al. 2000, but see also 
Ostriker et al. 2003). 

Determining the global matter density from cluster 
mass-to-light ratios therefore requires knowledge of the de- 
pendence of mass-to-light ratios on environment. Bahcall 
et al. (1995) show that mass-to-light ratios increase with 
scale from galaxies to groups to clusters. Ellipticals have 
larger overall values of M/Lb than spirals, presumably a 



result of younger, bluer stellar populations in spirals. At 
the scale of cluster virial radii, mass-to-light ratios appear 
to reach a maximum value. Some estimates of the mass-to- 
light ratio on very large scales (>10/i~"'^Mpc) are available 
(see references in Bahcall et al. 1995), but the systematic 
uncertainties are large. 

There are few estimates of mass-to-light ratios on scales 
between cluster virial radii and scales of lOft^^Mpc (Eisen- 
stein et al. 1997; SmaU et al. 1998; Kaiser et al. 2004; Rines 
et al. 2000, 2001a; Biviano & Girardi 2003; Katgert et al. 
2004; Kneib et al. 2003). On these scales, many galax- 
ies near clusters are bound to the cluster but not yet in 
equilibrium (Gunn & Gott 1972). These cluster infall re- 
gions have received relatively little scrutiny because they 
are mildly nonlinear, making their properties very diffi- 
cult to predict analytically. However, these scales are ex- 
actly the ones in which galaxy properties change dramati- 
cally (EUingson et al. 2001; Lewis et al. 2002; Gomez et al. 
2003; Treu et al. 2003; Balogh et al. 2004, and references 
therein). Variations in the mass-to- light ratio with envi- 
ronment could have important physical implications; they 
could be produced either by a varying dark matter fraction 
or by variations in the efficiency of star formation with en- 
vironment. In blue light, however, higher star formation 
rates in field galaxies compared to cluster galaxies could 
produce lower mass-to-light ratios outside cluster cores re- 
sulting only from the different contributions of young and 
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old stars to the total luminosity (Bahcall et al. 2000). 

Because clusters are not in equilibrium outside the virial 
radius, neither X-ray observations nor Jeans analysis pro- 
vide secure mass determinations at these large radii. There 
arc now two methods of approaching this problem: weak 
gravitational lensing (Kaiser et al. 2004) and kinematics 
of the infall region (Diaferio & Geller 1997; Diaferio 1999, 
hereafter DG97 and D99). Kaiser et al. (2004) analyzed 
the weak lensing signal from a supercluster at z « 0.4; 
the mass-to-light ratio (M/Lb=280 ±40 h for early-type 
galaxy light) is constant on scales up to 6 /i~^Mpc. Wil- 
son et al. (2001) finds similar results for weak lensing in 
blank fields; Gray et al. (2002) obtain similar results for a 
different supercluster. Recently, Kneib ct al. (2003) used 
weak lensing to estimate the mass profile of CL0024+1654 
to a radius of 3.25ft~^Mpc. Kneib et al. (2003) conclude 
that the mass-to-light ratio is roughly constant on these 
scales. 

Galaxies in cluster infall regions produce sharp features 
in redshift surveys (Kent & Gimn 1982; Shectman 1982; dc 
Lapparent et al. 1986; Kaiser 1987; Ostriker et al. 1988; 
Regos & Geller 1989). Early investigations of this infall 
pattern focused on its use as a direct indicator of the 
global matter density Unfortunately, random mo- 

tions caused by galaxy-galaxy interactions and substruc- 
ture within the infall region smear out this cosmological 
signal (DG97, Vedel & Hartwick 1998). Instead of sharp 
peaks in redshift space, infall regions around real clusters 
typically display a well-defined envelope in redshift space 
which is significantly denser than the surrounding environ- 
ment (Rines et al. 2003, hereafter Paper I, and references 
therein). 

DG97 analyzed the dynamics of infall regions with nu- 
merical simulations and found that in the outskirts of clus- 
ters, random motions due to substructure and non-radial 
motions make a substantial contribution to the amplitude 
of the caustics which delineate the infall regions (see also 
Vedel & Hartwick 1998, and references therein). DG97 
showed that the amplitude of the caustics is a measure 
of the escape velocity from the cluster; identification of 
the caustics therefore allows a determination of the mass 
profile of the cluster on scales < 10/i~^Mpc. 

DG97 and D99 show that nonparametric measurements 
of caustics yield cluster mass profiles accurate to ~50% 
on scales of up to 10 hr^ Mpc. This method assumes 
only that galaxies trace the velocity field. Indeed, simula- 
tions suggest that little or no velocity bias exists on linear 
and mildly non-linear scales (Kauffmann et al. 1999a,b). 
Geller et al. (1999, hereafter GDK), applied the kinematic 
method of D99 to the infall region of the Coma cluster. 
GDK reproduced the X-ray derived mass profile and ex- 
tended direct determinations of the mass profile to a ra- 
dius of 10 /i~^Mpc. The caustic method has also been 
applied to the Shaplcy Supercluster (Reisenegger et al. 
2000), A576 (Rines ct al. 2000, hereafter ROO), AWM7 
(Koranyi & Geller 2000), the Fornax cluster (Drinkwa- 
ter et al. 2001), A1644 (Tustin et al. 2001), A2199 (Rines 
et al. 2002), and six other nearby clusters (Paper I). Bi- 
viano & Girardi (2003) applied the caustic technique to an 
ensemble cluster created by stacking redshifts around 43 
clusters from the 2dF Galaxy Redshift Survey. ROO found 
an enclosed mass-to- light ratio of M/Lr 300/i within 



4 /i^iMpc of A576. Rines et al. (2001a) used 2MASS pho- 
tometry and the mass profile from GDK to compute the 
mass-to-light profile of Coma in the K-band. They found 
a roughly flat profile with a possible decrease in M/Lk 
with radius by no more than a factor of 3. Biviano & 
Girardi (2003) find a decreasing ratio of mass density to 
total galaxy number density. For early-type galaxies only, 
the number density profile is consistent with a constant 
mass-to-light (actually mass-to-number) ratio. 

Here, wc! calculate the infrared mass-to-light profile 
within the turnaround radius for the CAIRNS clusters 
(Paper I), a sample of nine nearby rich, X-ray luminous 
clusters. We use photometry from the Two Micron All 
Sky Survey (2MASS, Skrutskie et al. 1997) and add sev- 
eral new redshifts to obtain complete or nearly complete 
surveys of galaxies up to 1-2 magnitudes fainter than 
(as determined by Kochanek et al. 2001; Cole et al. 2001, 
hereafter KOI and COl). Infrared light is a better tracer of 
stellar mass than optical light (Gavazzi et al. 1996; Zibetti 
et al. 2002); it is relatively insensitive to dust extinction 
and recent star formation. Despite these advantages, there 
are very few measurements of infrared mass-to- light ratios 
in clusters (Tustin et al. 2001; Rines et al. 2001a; Lin et al. 
2003). 

Mass-to-light ratios within virial regions (where the 

masses arc more accurate than in the infall regions) pro- 
vide interesting constraints on the distribution of dark 
matter and stellar mass (see also Lin et al. 2003, here- 
after LOS). The virial masses in our sample span an order 
of magnitude in mass. More massive clusters have larger 
mass-to-light ratios. 

Cluster virial regions also provide potentially important 
constraints on the halo occupation distribution (e.g., Pea- 
cock & Smith 2000; Berfind & Weinberg 2002; Berfind 
et al. 2003, and references therein), the number of galax- 
ies in a halo of a given mass (see Cooray & Shcth 2002, 
for a recent review). The motions of galaxies and hot 
gas yield estimates of the dynamical mass independent 
of the number of galaxies (provided enough galaxies are 
present to yield a virial mass). Our mass profiles in Paper 
I arc among the first to extend significantly beyond r2oo- 
Thus, they should provide accurate estimates of r^oo- Also, 
the recent release of 2MASS allows us to count galaxies 
based on their near- infrared light, which is close to select- 
ing galaxies by stellar mass. Thus, both the masses and 
galaxy numbers are better defined than the few previous 
direct estimates of the halo occupation function (Peacock 
& Smith 2000; Marinoni & Hudson 2002; Lin et al. 2004). 

We describe the cluster sample, the near-infrared pho- 
tometry, and the spectroscopic observations in § 2. We 
discuss the galaxy properties (luminosity functions and 
broadband colors) within and outside the virial radius 
and compare both populations to field galaxies in §4. We 
calculate the number density and luminosity density pro- 
files in §5 and compare them to simple theoretical mod- 
els. We compute radial profiles of the mass-to-light ratio 
in §5. In §6 we constrain the halo occupation distribu- 
tion for the CAIRNS clusters and explore the dependence 
of mass-to-light ratios on halo mass. We discuss possible 
systematic uncertainties and the implications of our re- 
sults in §7 and conclude in §8. We assume a cosmology of 
Hq = 100 h km s~^,Ctm = 0.3, Oa = 0.7 except as noted 
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in §7. 

2. OBSERVATIONS 

2.1. The CAIRNS Cluster Sample 

We selected the CAIRNS parent sample from all nearby 
{cz < 15,000 km s"^), AbcU richness class R>1 (AbeU 
et al. 1989), X-ray luminous {Lx > 2.5 x 10*^h-'^erg 
s~^) galaxy clusters with declination 6 > —15°. Using X- 
ray data from the X-ray Brightest Abell Clusters catalog 
(EbeHng et al. 1996), the parent cluster sample contains 
14 systems. Wc selected a representative sample of 8 of 
these 14 clusters (Table 1). The cluster properties listed 
in Table 1 are from Paper I. The 6 clusters meeting the 
selection criteria but not targeted in CAIRNS are: A193, 
A426, A2G63, A2107, A2147, and A2657. The 8 CAIRNS 
clusters span a variety of morphologies, from isolated clus- 
ters (A496, A2199) to major mergers (A168, A1367). 

The redshift limit is set by the small aperture of the 1.5- 
m Tillinghast telescope used for the vast majority of our 
spectroscopic observations. The richness minimum guar- 
antees that the systems contain sufficiently large numbers 
of galaxies to sample the velocity distribution. The X-ray 
luminosity minimum guarantees that the systems are real 
clusters and not superpositions of galaxy groups (cf. the 
discussion of A2197 in Rines et al. 200ia, 2002). Three 
additional clusters with smaller X-ray luminosities (A147, 
A194 and A2197) serendipitously lie in the survey regions 
of A168 and A2199. A147 and A2197 lie at nearly identical 
redshifts to A168 and A2199; their dynamics are proba- 
bly dominated by the more massive cluster (Rines et al. 
2002). A194, however, is cleanly separated from A168 and 
we therefore analyze it as a ninth system. The inclusion of 
A194 extends the parameter space covered by the CAIRNS 
sample. The X-ray temperature of A194 listed in (Ebeling 
et al. 1996) is an extrapolation of the Lx — Tx relation; in 
Table 1 we therefore list the direct temperature estimate of 
Fukazawa et al. (1998) from ASCA data. Fukazawa et al. 
(1998) lists X-ray temperatures for 6 of the 8 CAIRNS 
clusters which agree with those listed in Ebeling et al. 
(1996). 

In Paper I, we applied a hierarchical clustering analysis 
(described in D99) to the redshift catalogs to determine 
the central coordinates and redshift of the largest system 
of galaxies in each cluster. Table 2 lists these hierarchi- 
cal centers and their projected separations from the X-ray 
peaks. We adopt these hierarchical centers as the cluster 
centers. 

2.2. 2MASS Photometry 

2MASS is an all-sky survey with uniform, complete pho- 
tometry (Nikolaev et al. 2000) in three infrared bands (J, 
H, and Kg, a modified version of the K filter truncated at 
longer wavelengths). We use photometry from the final 
extended source catalog (XSC, Jarrett et al. 2000). The 
2MASS XSC computes magnitudes in the Ks-hand using 
several different methods, including aperture magnitudes 
(using a circular aperture with radius 7"), isophotal mag- 
nitudes which include light within the elliptical isophote 
corresponding to /iif^=20 mag/arcsec^, Kron magnitudes, 
and extrapolated "total" magnitudes (Jarrett et al. 2000). 
The sky coverage of the catalog is complete except for 
small regions around bright stars. 



The 2MASS isophotal magnitudes omit ~15% of the to- 
tal flux of individual galaxies (KOI). COl compare 2MASS 
photometry from the Second Incremental Data Release 
(2IDR) with deeper infrared photometry from Loveday 
(2000). They find that Kron magnitudes are slightly 
fainter than the total magnitudes in deeper surveys (see 
also Andreon 2002a) and that 2MASS extrapolated total 
magnitudes are slightly brighter than Kron (roughly total) 
magnitudes from the deeper survey. 2MASS is a relatively 
shallow survey and thus likely misses many low surface 
brightness galaxies Andreon (2002a); Bell et al. (2003). 
In this work we focus on bright galaxies (which typically 
have high surface brightness) so this bias is less important 
than, e.g., estimates of the luminosity density or stellar 
mass density. 

Except where stated otherwise, we use the Ks-hand sur- 
vey extrapolated "total" magnitudes. Galactic extinction 
is usually negligible in the near-infrared. Wc correct for 
Galactic extinction by using the value in the center of the 
cluster. We make K corrections and evolutionary correc- 
tions of <0.15 magnitudes based on Poggianti (1997). Be- 
cause these corrections are small and not strongly depen- 
dent on the galaxy model at the redshifts of the CAIRNS 
clusters, we apply a uniform correction for all galaxies in a 
given cluster interpolated from the model Elliptical SED 
with solar metallicity and a star-formation e-folding time 
of 1 Gyr. 

Wc reprocess two galaxies in A576 and two galaxies near 
A2199 using the methodology of the 2MASS Large Galaxy 
Atlas (Jarrett et al. 2003). The galaxies in A576 (CGCG 
261-056 NEDOl and CGCG 261-056 NED 02) are bright 
ellipticals near the cluster center and also close to a bright 
star. One of the galaxies near A2199, UGC 10459, is an 
extremely flat edge-on disk galaxy. The other, NGC 6175, 
shows two nuclei aligned NW-SE. The SE component is 
brighter in Kg band. 

2.3. Spectroscopy 

The 2MASS photometry allows selection of complete, 

near- infrared-selected samples extending ^1-2 magnitudes 
fainter than the — —23.77 + 51og h determined for 

the field galaxy luminosity function in 2MASS extrapo- 
lated magnitudes (COl). We define i^r^-selected samples 
according to these magnitude limits within the smaller of 
the turnaround radius rt (the radius within which the av- 
erage density is 'S.bpc) or the limiting radius rmax of the 
caustic pattern (our membership criterion) in each cluster 
(see Paper I). Table 4 lists these radii and the apparent and 
absolute magnitude limits of these catalogs for the 9 clus- 
ters. Our redshift catalogs are 99.7% complete for cluster 
galaxy candidates brighter than Mk^ = — 23 -I- 51og/i and 
97.6% complete for candidates brighter than + 1. 

Most of the galaxies in these samples have redshifts 
in the redshift catalogs from the CAIRNS project (Pa- 
per I). Between 2002 June and 2003 September, we col- 
lected new redshifts for 515 galaxies with the FAST spec- 
trograph (Fabricant et al. 1998) on the 1.5-m Tilling- 
hast telescope of the Fred Lawrence Whipple Observatory 
(FLWO). FAST is a high throughput, long slit spectro- 
graph with a thinned, backside illuminated, antireflection 
coated CCD detector. The slit length is 180"; our observa- 
tions used a slit width of 3" and a 300 lines mm~^ grating. 
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Table 1 
CAIRNS Parent Population 



Cluster 


X-ray Coordinates 






Lx 


Tx 


Richness 




RA (J2000) 


DEC (J2000) 


km s~^ 


km 


10*3/1-2 ergs s-i 


keV 




A119 


00 56 12.9 


-01 14 06 


13268 


698^31 


8.1 


5.1 


1 


A168 


01 15 08.8 


+00 21 14 


13395 


579l|« 


2.7 


2.6 


2 


A496 


04 33 35.2 


-13 14 45 


9900 


701 +35 


8.9 


4.7 


1 


A539 


05 16 32.1 


-1-06 26 31 


8717 


734lij2 


2.7 


3.0 


1 


A576 


07 21 31.6 


-1-55 45 50 


11510 


1009+3g 


3.5 


3.7 


1 


A1367 


11 44 36.2 


-1-19 46 19 


6495 


782156 


4.1 


3.5 


2 


Coma 


12 59 31.9 


+27 54 10 


6973 


1042l|3 


18.0 


8.0 


2 


A2199 


16 28 39.5 


-1-39 33 00 


9101 


79611 


9.1 


4.7 


2 


A194 


01 25 50.4 


-01 21 54 


5341 


4951- 


0.4 


2.6 






Table 2 

CAIRNS Hierarchical Centers 



Cluster 


Hierarchical Center 
RA (J2000) DEC (J2000) 


km s"^ 


CZCMB 
km s^"'^ 


AR 


A119 


00 56 10.1 


-01 15 20 


13278 


12948 


56 


A168 


01 15 00.7 


-1-00 15 31 


13493 


13176 


239 


A496 


04 33 38.6 


-13 15 47 


9831 


9786 


24 


A539 


05 16 37.0 


-1-06 26 57 


8648 


8650 


33 


A576 


07 21 32.0 


-1-55 45 21 


11487 


11561 


16 


A1367 


11 44 49.1 


-1-19 46 03 


6509 


6837 


61 


Coma 


13 00 00.7 


-1-27 56 51 


7096 


7365 


153 


A2199 


16 28 47.0 


-1-39 30 22 


9156 


9181 


86 


A194 


01 25 48.0 


-01 21 34 


5317 


5011 


11 



This setup yields spectral resolution of 6-8 A and covers 
the wavelength range 3600-7200 A. We obtain redshifts 
by cross-correlation with spectral templates of emission- 
dominated and absorption-dominated galaxy spectra cre- 
ated from FAST observations (Kurtz & Mink 1998). The 
typical uncertainty in the redshifts is 30 km s~^. Table ?? 
lists the new redshifts. The additional redshifts make no 
significant difference to the locations of the caustics or to 
the resulting mass profiles. We thus use the caustics and 
mass profiles from Paper I. 

An important difference between the FAST spectra col- 
lected for CAIRNS and those collected for other, larger 
redshift surveys (Colless et al. 2001; Stoughton et al. 2002) 
is that CAIRNS suffers no incompleteness due to fiber 
placement constraints. 

We calculate the maximum fraction fnoz of light missing 
from our catalogs if we assume that all galaxies without 
redshifts and brighter than the magnitude limit are clus- 
ter members (Table 4). In other words, we evaluate the 
potential observational bias which results if every galaxy 
without a redshift were a cluster member. For this extreme 
case, the total luminosity within rmax is underestimated 
by the fraction fnoz < 0.10 for all clusters. The new red- 
shifts in Table ?? contribute significantly to the complete- 
ness of these catalogs. Because the galaxies without red- 
shifts are almost entirely faint galaxies at large distances 
from the cluster center, fnoz is a very conservative upper 
limit on the fraction of light missing within the complete- 
ness limits (the surface number density of member galax- 
ies decreases with radius and the fraction of background 



galaxies increases with apparent magnitude). 

Assuming that the luminosity function in clusters and 
infall regions is identical to that in the field (we test this 
assumption in §3.1), we calculate the fraction of total 
galaxy light contained in galaxies brighter than our com- 
pleteness limits. This fraction is greater than 60% for all 
clusters. From repeated measurements, apparent magni- 
tudes in the 2MASS XSC have an uncertainty of ~ 0.14 
magnitude at Kg = 13.4 — 13.5; the galaxy catalogs prob- 
ably suffer incompleteness fainter than Kg « 14. Thus, 
the 2MASS XSC provides accurate magnitudes for galax- 
ies within our completeness limits, but it is difficult to 
use 2MASS galaxy counts at much fainter magnitudes to 
estimate the contribution of fainter galaxies to the total 
cluster/infall region light. Note that the field luminosity 
function of COl that we adopt here has a steeper faint-end 
slope than the luminosity function calculated from Kron 
magnitudes. If we adopt the Kron magnitude faint-end 
slope of COl, /l increases by 7-15% (the best sampled 
clusters have the smallest changes). We discuss this issue 
further in §3.1 and §7.2. 

Figure 1 shows the redshift completeness as a function of 
apparent and absolute magnitude {Ks extrapolated mag- 
nitude) along with the total number of galaxies, the num- 
ber with redshifts, and the number of members versus 
magnitude. Note that, as in Paper I, we order the clusters 
by decreasing X-ray temperature from left to right and 
from top to bottom in this and all similar later figures. 
The catalogs are complete for cluster galaxy candidates 
brighter than Mk^ = — 23 except for five candidates in 
the outskirts of A539 which lie at high Galactic extinc- 
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tion. It is not clear whether these objects are galaxies 
or extended Galactic infrared sources. The brightest of 
these sources, IRAS 05155+0707, is an embedded Class 
1 protostar and likely the source of Herbig-Haro objects 
HH114 and HH115 (Reipurth et al. 1997). We exclude 
IRAS 05155+0707 from the photometric catalog and the 
calculation of fnoz in Table 4. 

Figure 1 also shows constraints on the luminosity func- 
tions in the clusters. The sets of dash-dotted lines show the 
limits from assuming that (1) all galaxies without redsliifts 
are members or (2) none are. We discuss the luminosity 
functions in more detail in §4, but we note here that the 
faint-end slope of the luminosity function in infall regions 
is poorly constrained without deep, complete spectroscopy. 

3. PROPERTIES OF GALAXIES INSIDE AND OUTSIDE THE 
VIRIAL REGION 

Galaxy properties such as morphology and star forma- 
tion rate are strongly correlated with their local and global 
environments (e.g., EUingson et al. 2001; Lewis et al. 2002; 
Gomez et al. 2003; Trcu et al. 2003; Balogh et al. 2004, 
and references therein). Differences in galaxy properties 
with environment may lead to apparent changes in the 
observed mass-to-light ratio even if the ratio of dark mat- 
ter to stellar mass remains constant (e.g. Bahcall et al. 
2000). The CAIRNS 2MASS selected galaxies provide a 
well-defined population with which to investigate these en- 
vironmental effects. The environments considered range 
from cluster centers with densities ^lOOOpc to the edges of 
infall regions with densities ^3pc at the turnaround radius 
rt- These environments are all denser than the universal 
average density ^mpc, but they cover the density range 
where galaxy morphologies, optical colors, and star for- 
mation rates change dramatically (EUingson et al. 2001; 
Lewis et al. 2002; Gomez et al.' 2003; Treu et al. 2003; 
Balogh et al. 2004). We investigate the near-infrared pho- 
tometric properties of galaxies inside and outside the virial 
regions of the CAIRNS clusters and compare them to field 
galaxies. 

3.1. Luminosity Functions 

Many investigators have sought to determine the en- 
vironmental dependence of the luminosity function (e.g., 
Balogh et al. 2001; Beijersbergen et al. 2002; De Propris 
et al. 2003, and references therein). Using the 2dF Galaxy 
Redshift Survey, De Propris et al. (2003) fit their cluster 
data to the Schechter (1976) luminosity function (LF), 

N{M) OC 100-4(«+l)(M*-M)g-loO-("*-") 

and find that the cluster LF in the bj band has a brighter 
characteristic magnitude M* and steeper faint-end slope 
a than the field LF. Although differences between cluster 
and field luminosity functions exist at other wavelengths 
(e.g., Trentham 1998a,b; Mobasher et al. 2003; Sabatini 
et al. 2003), the cluster LF in the Kg band is quite similar 
to the field LF (e.g., Mobasher & Trentham 1998; de Pro- 
pris et al. 1998; Andreon & Pello 2000; Tustin et al. 2001; 
Balogh et al. 2001), perhaps indicating a universal stellar 
mass function (Andreon 2004). Balogh et al. (2001) com- 
bine data from 2MASS and the Las Campanas Redshift 
Survey and find that the cluster LF in the J band has a 
brighter characteristic magnitude and a steeper faint-end 



slope than the field LF; similar differences are seen at Ks 
band but the parameters differ by <3-(j. Andreon (2004) 
finds that the cluster and field LFs are indistinguishable 
at red wavelengths in the optical (see also Christlein & 
Zabludoff 2003) , suggesting that much of the difference at 
bluer wavelengths is due to star formation. 

Figure 2 shows the near-infrared luminosity functions of 
each of the CAIRNS clusters including all galaxies within 
the infall regions. We use the caustics from Paper I to 
define membership. In magnitude bins without complete 
redshifts, we compute a completeness correction by assum- 
ing that the membership fraction of galaxies without red- 
shifts is the same as the membership fraction of galaxies 
with redshifts. Galaxies without redshifts tend to be at 
larger projected clustrocentric distances than those with 
redshifts. One might thus expect that these galaxies are 
more likely to be non-members because the ratio of clus- 
ter members to background galaxies decreases with radius. 
Counteracting this effect, galaxies without redshifts tend 
to have lower surface brightnesses than those with red- 
shifts (because of observational bias towards higher sur- 
face brightness galaxies); because of the correlation be- 
tween absolute magnitude and surface brightness, galaxies 
of a given apparent magnitude with lower surface bright- 
nesses should be intrinsically fainter and are thus more 
likely to be cluster members (Conselice et al. 2002; Ko- 
ranyi & Geller 2000). 

We count the number of bright galaxies (those with 
Mk^ < + 1) in each cluster and use this number 

to calculate relative normalizations for each cluster. Fig- 
ure 2 shows the Schechter LF for field galaxies from COl 
scaled by this number of bright galaxies with an arbitrary 
overall normalization. 

We compute the luminosity functions separately for the 
virial regions and the infall regions taking i?200 {fs is the 
radius within which the enclosed mass density is S times 
the critical density, Rs — rs is the projected radius) as 
the dividing radius. Some galaxies projected inside R200 
lie outside r200j but no galaxies projected outside i?200 
lie inside r2oo; thus the luminosity functions inside i?200 
will be contaminated by galaxies outside the virial region. 
Figure 3 shows the luminosity functions within i?2oo ; Fig- 
ure 4 shows the luminosity functions outside this radius. 
In each panel, we plot the best-fit Schechter (1976) lumi- 
nosity function for field galaxies from COl scaled by the 
number of bright galaxies with an arbitrary overall nor- 
malization. Figure 5 shows the combined CAIRNS LFs 
inside and outside R200 as well as the total LF. The LFs 
in the virial regions and infall regions are very similar. 

At the bright end, the LFs in both the virial regions 
and the infall regions are poorly fit by a Schechter func- 
tion (Figure 5); the observed LFs contain more galax- 
ies brighter than M/f, = —25 and fewer galaxies at 
—25 < Mk^ < —24 than a Schechter function which fits 
the faint-end slope. This difference may result from the 
existence and evolution of cD galaxies (e.g., Schombert 
1988; Tonry 1987) present only in cluster environments. 
Figure 6 shows the ratio of the LF outside i?200 to that 
inside i?2oo • The infall region LF contains fewer extremely 
bright galaxies {Mk^ < —25) than the virial region LF, 
but there is very little difference within the Poissonian un- 
certainties. Also, it is worth noting that extremely bright 
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galaxies arc present in six of the nine infall regions (Coma, 
A119, A2199, A576, A168, and A194), demonstrating that 
these bright galaxies do not reside exclusively in cluster 
centers. Many of these galaxies likely occupy the centers 
of galaxy groups in the infall regions (Rines et al. 2001b). 
A test shows that the LF ratios for all bright galaxies 
(Afifa < —23.7) are consistent with a constant value at the 
95% confidence level. 

At magnitudes fainter than the completeness limit, the 
LF in the infall region (excluding the virial region) consis- 
tently exceeds that inside the virial region, suggesting that 
the faint-end slope might be steeper in the infall region. 
The uncertainties in Figure 6 are Poissonian. Because the 
correction for galaxies without redshifts may be biased, 
these uncertainties may be significantly underestimated. 
A deeper complete spectroscopic survey of the infall re- 
gions is necessary to determine the reality of effects at 
these faint luminosities. 

We calculate the best-fit luminosity function of the 
Schechter (1976) form for Mk, < -22.1 -|- 5\ogh for ah 
the clusters combined. This limit corresponds to the 
2MASS completeness limit of Ks = 13.5 for the most dis- 
tant CAIRNS clusters. We fit the LF for galaxies within 
.R2001 outside i?200, and all galaxies combined. We do not 
account for measurement uncertainties in the fits. Ta- 
ble 5 lists the best-fit parameters (from minimizing x^) 
as well as determinations of the field luminosity function 
(KOljCOl). The uncertainties are 68% confidence limits 
for two interesting parameters. We list two different es- 
timates from COl, one using extrapolated magnitudes (as 
used here) and one using 2IDR Kron magnitudes converted 
to 'total' magnitudes by subtracting -0.20 mag (see COl). 
The LF parameters dififer by 2-3<t from the field values, 
and agree well with previous determinations (Balogh et al. 
2001, L03). However, the fits to the CAIRNS LPs are not 
very good; the probability of obtaining a larger value of 
from a sample drawn from the Schechter LF is <0.7% for 
the total LF. The best-fit characteristic magnitude of the 
virial region LF is brighter than the field LF, similar both 
in sign and magnitude to the difference found by Balogh 
et al. (2001); the faint-end slope of the CAIRNS virial re- 
gions is slightly steeper than the field values. The LFs in 
the infall regions are intermediate between the field LFs 
and the virial region LFs. 

We repeat the fits using the completeness limit of the 
redshift catalogs Mk, < —22.7 and obtain consistent pa- 
rameters with larger uncertainties due to the weaker con- 
straints on the faint ends of the LFs. We experimented 
with different cuts in absolute magnitude both at the 
bright end (excluding cD-like galaxies that could skew the 
LF parameters) and the faint end. The best-fit LF parame- 
ters are fairly sensitive to the limiting magnitude adopted, 
perhaps because the cluster LF is not well described by a 
Schechter function. However, these parameters are gen- 
erally within the 2-a range of the values listed in Table 
5. 

It is interesting that the characteristic magnitude of the 
CAIRNS virial region LF agrees well with that of the clus- 
ter LF constructed by LOS without spectroscopy. This 
agreement suggests that statistical background subtrac- 
tion produces little bias in the resulting LF parameters. 
Both L03 and CAIRNS use 2MASS photometry which pro- 



vides only a limited probe of the faint-end slope. It would 
be instructive to compare the LFs of individual clusters 
constructed with spectroscopic membership to those con- 
structed with statistical background subtraction. A de- 
tailed comparison is outside the scope of the present work, 
but in §5 we will show that LFs constructed with statis- 
tical background subtraction (L03) yield mass-to-light ra- 
tios consistent with our results for clusters with complete 
spectroscopy. 

The best-fit LF parameters significantly affect the esti- 
mates of the fraction of light Jl contained in faint galaxies 
(see Table 4). However, for fixed LF parameters, the ra- 
tio of the maximum to the minimum values of Jl for the 
clusters varies by less than 10%; thus, the relative values 
of /l are robust. Because the CAIRNS LF parameters are 
consistent with the field LF but have larger uncertainties, 
we continue to use the field LF to estimate the fraction 
of light contributed by faint galaxies. Note that the field 
LF we adopt (COl extrapolated magnitudes) has both a 
brighter characteristic magnitude and a steeper faint-end 
slope than the LF of COl from Kron magnitudes. 

We repeat this analysis in the J band, which extends 
deeper in 2MASS and thus has smaller statistical uncer- 
tainties. Figure 7 shows the J band LF for all galaxies 
within rt, Figure 8 shows the luminosity functions within 
i?200, and Figure 9 shows the luminosity functions outside 
this radius. We combine the LFs to produce an average 
cluster LF in Figure 10. We scale the LF inside and out- 
side i?200 to have the same normalization at M} for field 
galaxies. As in Kg band, the cluster LF has a very sim- 
ilar shape to the field LF except for an excess of bright 
galaxies. We repeat the non-parametric test of comput- 
ing the LF ratios (Figure 11). Table 5 lists the best-fit 
Schechter function parameters. These parameters differ 
by no more than 3-cr from the field values determined by 
COl. The characteristic magnitude Mj for cluster virial 
regions is brighter than the field value by about 0.5 mag- 
nitudes, consistent with the results of Balogh et al. (2001). 
There is remarkably little difference between the two LFs 
across the entire range of magnitudes, although at faint 
magnitudes there is room for significant differences which 
could be explored with deep, complete spectroscopy. 

To summarize, we see marginal evidence for differences 
between the cluster LF and the field LF. The cluster LF 
is slightly brighter and has a steeper faint-end slope than 
the field LF. We obtain similar results in both J and Ks 
bands. Our data sample only giant galaxies, so signifi- 
cant differences may exist in the cluster and field LFs in 
the dwarf galaxy regime. For the purposes of comput- 
ing mass-to-light ratios, the systematic uncertainty intro- 
duced by possible differences in the cluster and field LFs is 
< 15%. Note that, as expected, the LF in the infall region 
is intermediate between the field LF and the cluster LF. 

3.2. Luminosity Segregation 

Dynamical friction could lead to luminosity segregation 
in galaxy clusters. Some investigators have claimed evi- 
dence for luminosity segregation in compilations of cluster 
data (e.g. Adami et al. 1998; Andreon 2002b, and refer- 
ences therein). Figure 12 shows the distribution of abso- 
lute magnitude versus (projected) distance Rp from the 
cluster center. If luminosity segregation were significant. 
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we would see more bright galaxies near cluster centers. 
The brightest cluster galaxy is typically very close to the 
cluster center, consistent with a bright central cD galaxy 
increasing in mass through accretion of smaller galaxies. 
However, there arc also many comparably bright galaxies 
in the outskirts of the clusters. In A2199, many of the 
extremely bright galaxies outside the virial region are at 
the centers of infalling groups (Rines et al. 2001b, 2002). 
There is little evidence for luminosity segregation in the 
CAIRNS clusters, consistent with earlier results for A576 
(Rines et al. 2000). This result is not surprising given 
the similarity of the LFs inside and outside i?2oo (Figure 
5). Again, note that the CAIRNS samples do not extend 
into the dwarf galaxy regime, where luminosity segregation 
might be present (Andreon 2002b). 

3.3. Broadband Galaxy Colors 

Star formation rates depend on environment (e.g., 
EUingson et al. 2001; Lewis ct al. 2002; Gomez et al. 
2003; Treu et al. 2003; Balogh et al. 2004, and refer- 
ences therein). Because stellar populations in field galax- 
ies are on average younger than those in cluster galax- 
ies, more blue light is emitted per unit mass in field-like 
environments than in cluster environments. As a conse- 
quence, mass-to-blue-light profiles might decrease with ra- 
dius (Bahcall et al. 2000) even if the ratio of gravitational 
mass to stellar mass were constant. 

Because young stars are both hotter and bluer than 
older stars, the difference in stellar mass-to-light ratios 
decreases toward longer wavelengths (see the synthesized 
stellar population models of Bruzual & Chariot 2003). For 
example, studies of near-infrared mass-to-light ratios in 
galaxies suggest that the mass-to-light ratio at these wave- 
lengths is insensitive to the current star formation rate 
in either disk galaxies (Gavazzi et al. 1996) or early-type 
galaxies (Zibetti et al. 2002). Unfortunately, the color dif- 
ferences between J and K bands arc not very large because 
these wavelengths primarily trace Population II stars ( Jar- 
rett et al. 2003), making this effect difficult to detect with 
2MASS data alone. 

For A576, our 9 square degrees of photometric CCD ob- 
servations in R band (Rines et al. 2000) allow a measure- 
ment oi R — Kg colors. Optical- infrared colors enable us to 
investigate stellar population effects. Although both the 
2MASS magnitudes used here and the R band magnitudes 
in Rines et al. (2000) are supposed to be close to total, a 
systematic difference in the magnitude definitions could 
introduce an artificial color gradient. We reprocess the R 
band images using SExtractor (Bertin & Arnouts 1996) 
to obtain aperture magnitudes within a circular aperture 
of radius 14" in R band and radius 15" in 2MASS. This 
slight mismatch in apertures produces a small bias towards 
redder colors, but clustrocentric gradients, if any, should 
still be evident. We calculate R— J and R— Kg colors for 
all of the galaxies in both catalogs. Figure 13 displays the 
R—Ks colors versus projected radius. There is no obvious 
radial gradient in either R — J ov R — Kg colors for bright 
galaxies. There may be a radial gradient in R — Kg col- 
ors for galaxies fainter than Mk^ — —22.77, but we lack 
complete spectroscopy at these magnitudes. In a photo- 
metric study of clusters using SDSS, Goto et al. (2004) 
find small but significant radial gradients in the fraction 



of blue galaxies with radius (the fraction increases with 
radius). The trends are weakest in the most nearby clus- 
ters which are the most similar to the CAIRNS clusters. 
Note that the CAIRNS catalogs are selected at Kg rather 
than r, which may account for the lack of a gradient in 
A576. Also, the trends may be weaker m. R — Kg colors 
than in, e.g., u — r colors, which are much more sensitive 
to the presence of young stars. A multiwavelength study 
of several clusters with spectroscopically determined mem- 
bership would clarify the importance of color gradients in 
clusters. 

We plot the R—Kg color versus Kg magnitude in Figure 
14. There is little evidence for a color-magnitude relation. 
Galaxies inside and outside i?2oo occupy the same parts 
of the diagram, indicating that there is no large difference 
in the two populations. The galaxies appear to have very 
similar stellar populations. Comparing the colors to the 
models of Bruzual & Chariot (2003) indicates metallici- 
ties greater than solar. The degeneracy between age and 
metallicity effects prevents further conclusions. 

3.4. Near-Infrared Colors and the Color- Magnitude 
Relation 

Significant variations in the stellar mass-to-light ratio 
might be indicated by radial gradients in J — Kg colors. 
Unfortunately, identifying such gradients is difficult be- 
cause variations in galaxy J — Kg colors are relatively small 
and because unlike optical colors, galaxies with the reddest 
J — Kg colors may contain younger stellar populations and 
are red as a result of emission from hot dust (Hunt et al. 
2002; Barton Gillespie et al. 2003). We observe no radial 
color gradients in J — Kg (sec Figure 15, which highlights 
the lack of trends in the outlying points). Figure 16 shows 
that the distributions oi J — Kg colors of bright galaxies 
inside and outside J?200 are extremely similar. There is 
a possible excess of galaxies in the red tail of the distri- 
bution in the sample outside i?2oo (Figure 15); some of 
these are edge-on disk galaxies while others are probably 
AGN (Jarrett 2000; Jarrett et al. 2003). Because of the 
morphology-density relation, we expect more disk galaxies 
in cluster outskirts. 

We plot J — Kg color (computed within the ellipti- 
cal isophote Kg=2Q mag arcsec"^) versus absolute mag- 
nitude Mk, (from the extrapolated Kg magnitude) of 
CAIRNS members in Figure 17. The most striking re- 
sult is that the outlying data points are galaxies both 
inside and outside i?200, which suggests that the stellar 
populations of galaxies in these regions are similar. There 
is tentative evidence for a color-magnitude relation (i.e., 
fainter galaxies arc bluer, see, e.g., Terlcvich ct al. 2001, 
and references therein) in the near-infrared, but the slope 
(w — O.Olmag mag^^) is much shallower than at optical 
wavelengths, e.g., the slope is — 0.14±0.01 in U—V versus V 
in Coma (Terlevich et al. 2001). We obtain a similar color- 
magnitude relation when the colors and magnitudes are 
determined from aperture photometry, suggesting that the 
relation does not result from systematic effects in 2MASS. 
The variations in the colors can be explained by variations 
in the metallicities of the stellar populations. The most re- 
cent stellar population models of Bruzual & Chariot (2003) 
indicate that 10 Gyr old stellar populations (formed in- 
stantaneously according to a Chabrier (2003) initial mass 
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function) with mctallicities [Fe/H] = —0.64 +0.56 
have rest-frame J — Kg = 0.75 — > 1.1. As at optical wave- 
lengths, there is degeneracy between age and metallicity 
effects; bluer colors result from either lower mctallicities or 
younger ages (Worthcy 1994). Accurate spectral informa- 
tion is required to break this degeneracy (e.g., Concannon 
et al. 2000). 

The scatter in the observed near-infrared color- 
magnitude relation is larger for fainter galaxies; the fainter 
galaxies have more varied stellar populations and/or larger 
uncertainties. Note that galaxies in clusters at low galactic 
latitude (A539 and A496) have larger scatter than those in 
clusters near the galactic poles (Coma and A1367). This 
observation suggests that a significant part of the scatter 
may result from uncertainties in Galactic extinction. A 
full accounting of the color-magnitude relation is beyond 
the scope of this paper. Instead, we note that the near- 
infrared properties of galaxies do not change dramatically 
with radius. This result implies that the stellar mass-to- 
light ratios do not change dramatically with radius; thus, 
measuring near-infrared mass-to-light ratios is a good ap- 
proximation to a measurement of the ratio of total mass 
to stellar mass. 

4. NEAR-INFRARED LUMINOSITY AND NUMBER DENSITY 

PROFILES 

4.1. Number Density Profiles 

Because our catalogs are essentially complete within 
their respective magnitude limits, we can count the num- 
ber of bright galaxies to compare cluster richness. We 
adopt Mk, = —22.77 -|- 51og/i as our limiting magnitude 
because all clusters are complete to this depth (Table 4). 
This limit is equivalent to + 1 for field galaxies. Ta- 
ble 6 lists the number of galaxies inside and outside i?2oo 
("outside i?20o" means the projected radius Rp satisfies 
''200 < Rp < Tmax)- In all clustcrs but A539, there are 
more cluster members outside i?200 than inside i?2oo- We 
suggested this result in Paper I but lacked the uniform 
photometry necessary to establish it. 

Figure 18 shows the surface number density profiles of 
the CAIRNS clusters. We choose radial bins spaced loga- 
rithmically by 0.25; the outermost bin contains the maxi- 
mum radius r^ax of the caustics. We fit the number den- 
sity profiles of the CAIRNS clusters to three simple ana- 
lytic models. The simplest model of a self-gravitating sys- 
tem is a singular isothermal sphere (SIS). The volume den- 
sity of the SIS decreases with radius according to p oc r~^; 
the projected number density of objects S decreases as 
S(iip) oc R-^. Navarro et al. (1997) and Hernquist (1990) 
propose two-parameter models based on CDM simulations 
of haloes. These density profiles are 

p(r)<x[^(l + ^)-«] (2) 

where a is a scale radius and a=2 for the NFW profile 
and a=3 for the Hernquist profile. At large radii, the 
NFW density profile decreases as r~^ and the density of 

the Hernquist model decreases as r^^ (implying a finite 
total mass). The NFW surface number density profile is 

^'<''' ,ln(OT-l) ''-^'''l 
where s = Rp/a is the projected radius in units of the scale 
radius, N{a) is the number of galaxies within the sphere 



of radius a, and 



We fit the parameter N{a) rather than the core density 

ria = 3N{a)/ {Aira^) because N{a) and a arc much less 
correlated than ria and a (Mahdavi et al. 1999). The Hern- 
quist surface density profile is 

^) = ^?^[(2 + ^')^(^) - 3] (5) 
7ra^(s^ — 1)'' 

where a is the scale radius and M is the total mass. Note 
that M(a) = M/4. We minimize and list the best-fit 
parameters ajv for the NFW and Hernquist models in Ta- 
ble 7. We perform the fits on all data points within the 
maximum radii listed in Table 4. We plot the surface num- 
ber density profiles and the best-fit NFW (solid lines) and 
Hernquist (dash-dotted lines) models in Figure 18. The 
SIS (dashed lines) is not normalized and is shown only for 
comparison. 

The best-fit scale radii un for both the NFW and Hern- 
quist models arc larger than the best-fit scale radii cm of 
the mass profiles in Paper I for all clusters except A539, 
where the NFW scale radius is the same. In individual 
clusters, and differ only at l-3cr significance. How- 
ever, a K-S test indicates that the distributions of cat and 
om are not drawn from the same population at the 99.6% 
confidence level for the NFW model and at 99.95% confi- 
dence for the Hernquist model. These differences suggest 
that mass is more concentrated than light. 

Two clusters, A168 and A1367, are poorly fit by Hern- 
quist and NFW profiles. They both could be fit by these 
profiles within -R200! but the surface number density out- 
side i?2oo exceeds the predicted profile. This result sug- 
gests that these chistcrs arc not isolated from surrounding 
large-scale structure and that we may be observing them 
at an early stage of their evolution. In fact, both these 
clusters contain major mergers. Furthermore, they are the 
only CAIRNS clusters currently undergoing major merg- 
ers. Excluding these clusters from the comparison of scale 
radii only reduces the significance of the K-S test to 99.5% 
for both models. Thus, the difference in distribution of 
scale radii is not solely a result of these merging systems. 
Similarly, A2199 has a large core component. Excluding 
the innermost bin slightly increases the best-fit value of 
UN and slightly decreases the best- fit value of N{a). 

Biviano & Girardi (2003) find similar results from a 
Jeans analysis of an ensemble cluster constructed from the 
2dFGRS: the ratio of mass density to galaxy (deprojected) 
number density decreases with radius. Similarly, Lin et al. 
(2004) find that the concentration of galaxies is smaller 
than the expected concentration of mass, i.e., the galaxies 
arc more extended than expected. From our results and 
the independent analyses of these other authors, we thus 
conclude that the difference is a real physical effect. 

4.2. Luminosity Profiles 

Because the cluster LF is not significantly different from 
the field LF, the estimates of the fraction of the total 
cluster light contained in galaxies brighter than the magni- 
tude limits in Table 4 (which assume the LF parameters of 
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the field LF) are justified. We estimate the total light by 
adding the luminosity in galaxies brighter than the magni- 
tude limits, then dividing by /l. We make no corrections 
for the small incompleteness in our spectroscopic catalogs 
ifnoz in Table 4). This omission could lead to slight under- 
estimates of the luminosity in the outskirts of the clusters. 
The photometric uncertainties in the luminosity profiles 
are <10%. Because we compute the luminosity profiles 
only from relatively bright galaxies, the uncertainties are 
dominated by counting statistics (Kochanck et al. 2003). 

We fit the luminosity density profiles of the CAIRNS 
clusters to the simple analytic models described in the 
previous section, replacing N{a) with L(a), the luminos- 
ity contained within the sphere of radius a^. Figure 19 
shows the surface luminosity density profiles and the best- 
fit NFW (solid lines) and Hernquist (dash-dotted lines) 
models. The scale radii oj, of the light distributions are 
close to flTv and larger than ajv/, again implying that the 
light in galaxies is more extended than the mass. A K-S 
test indicates that the distributions of az, and um are not 
drawn from the same population at the 98.1% confidence 
level for the NFW model and at 99.95% confidence for 
the Hernquist model. A K-S test detects no differences in 
the distributions of az, and ajv for either model. Again 
we conclude that the mass is more concentrated than the 
light. 

5. NEAR-INFRARED MASS-TO-LIGHT PROFILES 

We next compute M(< r)/LK^{< R) as a function of 
radius (in units of Mq/Lq) using the caustic mass profiles 
and the luminosity profiles from the previous section (solid 
lines in Figure 20 with uncertainties shown in shading). 
The luminosity profiles Lk^{< R) are projected in two 
dimensions; the mass profiles M(< r) are radial profiles. 
Thus, these mass-to-light profiles are the mass in spheres 
divided by the light in cylinders. Table 8 summarizes the 
mass-to-light ratios inside and outside r2oo calculated by 
dividing the caustic masses (in spheres) by the light pro- 
files (in cylinders). Without correcting for this geometric 
effect, the mean mass-to-light ratio inside r2oo (70 ± 7h) 
is a factor of 1.8 ± 0.3 larger than the mean mass-to-light 
ratio outside r2oo (38±6/i). The mean mass-to- light ratio 
inside the maximum radius probed, Vmax, is 53 ± 5h. 

One notable feature of Figure 20 is the variety in the 
shapes of the mass-to- light profiles in individual clusters. 
Some clusters have flat profiles while others have strongly 
peaked profiles. There is no obvious cause of these differ- 
ences (e.g., presence of a bright cD galaxy, presence of a 
major merger). Indeed, significant variation in the shapes 
of mass-to-light profiles obtained with the caustic tech- 
nique is expected from projection effects along different 
lines of sight (D99). 

Obviously, it is preferable to compute both the mass 
and the light in either spheres or cylinders but not one 
of each. Becaiise mass and light are never negative, 
LkA< r) < LkA< Rp) and M{< r) < M{< Rp). Thus, 
this geometric effect should artificially decrease the "ob- 
served" mass-to-light ratios in the centers of the clusters. 
To correct for this geometric effect, one must make as- 
sumptions about the shapes of the profiles and their sym- 
metries. We prefer to present the data with few manipu- 
lations. We thus project the mass profiles into cylinders 



rather than attempting to deproject the noisy luminosity 
profiles. In particular, we assume that the mass distribu- 
tion is well-described by one of the simple mass models. 
Because the Hernquist profile has a finite total mass, the 
best-fit Hernquist profiles have smaller densities than the 
best-fit NFW profiles at large radii. Thus, the projected 
Hernquist profile is more centrally concentrated than the 
NFW profile. These projection effects will thus be larger 
if the true profile is a Hernquist profile. We show these 
projected Mh{< Rp) / LkX"^ R) niass-to-light profiles as 
dash-dotted lines in Figure 20. As expected, these pro- 
files have larger mass-to-light ratios at small radii than 
the spheres-by-cylinders profiles. If the decreasing shapes 
of these profiles are correct, the deprojected mass-to-light 
profiles M(< r)/L{< r) should decrease slightly faster 
than in projection. 

We quantify the size of this effect for NFW mass pro- 
files. For NFW profiles with concentrations c = 5 — 20 
(in Paper I we measured c = 5 — 17 for the CAIRNS 
clusters), the projected mass within a cylinder of radius 
i?200 is a factor of 1.15-1.25 larger than the mass in a 
sphere of radius r2oo (the factor increases with decreas- 
ing concentration). Projection effects are less dramatic for 
cylindrical shells compared to spherical shells at large radii 
because some light /mass outside the spherical shell is pro- 
jected into the cylindrical shell; some light/mass within the 
spherical shell is projected into cylindrical shells at smaller 
radii. For NFW halos with concentrations c = 5 — 20, 
the projected mass in the cylindrical shell bounded by 
ii200 and Rmax is 1-5% greater than the mass in the 
spherical shell boimded by r2oo and r^nax- Thus, if the 
CAIRNS clusters are well-described by NFW profiles with 
c «5, the mass-to-light ratio inside the cylinder R2Q0 is 
larger by a factor of 1.15-1.25 than the measured quan- 
tity M(< r2m)/L{< i?20o)- Similarly, the mass-to-light 
ratio in the cylindrical shell bounded by i?,200 and Rmax is 
larger by a factor of 1.01-1.05 than the measured quantity 
) - M{< r2oo)]/[-^^(< R max 

)-L{< i?2oo)]. The 
difference in mass-to-light ratios inside and outside R200 
is therefore larger than calculated above; 

(M/L)(< j?20o) ^^^ Mi<r2oo)/L{<R2oo) 
{M/L){> R200) ^ ■ M(> r2oo)/i(> R2Q0) ~ ■ 

(6) 

The projection effects therefore aggravate the difference 
in mass-to-light ratios between cluster virial regions and 
their outskirts. 

The preceding calculation used the nonparametric mass 
profiles from Paper I. We repeat the calculations of Ta- 
ble 8 using the best-fit (parametric) NFW mass profiles 
from Paper I projected into cylinders. The mean mass- 
to-light ratio inside r2oo is M^fw / L\r2oo = 77 ± 7/i and 
the mean mass-to-light ratio outside r2oo is 48 ± 5/i. The 
ratio of these is 1.6±0.2, very similar to the ratio 1.8±0.3 
calculated above with no corrections for geometric projec- 
tion effects. We note here that the best-fit NFW param- 
eters do not vary significantly if the fits are restricted to 
r < 1.5 /i~^Mpc; at these radii, the mass profiles agree 
with X-ray and virial mass estimates (see Paper I). 

We now calculate the mass-to-light profile in individual 
shells. The enclosed mass-to-light profiles calculated above 
decrease with radius. Because these profiles are dominated 
by the mass-to-light ratio in the core, the mass-to-light ra- 
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tio in the outer shells must generally be smaller than the 
enclosed mass-to- light ratio at that radius. The uneer- 
tainties in individual shells are sufficiently large that we 
must bin several shells to obtain a significant signal. These 
mass-to-light ratios are the mass in spherical shells divided 
by the light in cylindrical shells. As noted above, for spher- 
ically symmetric NFW models, the projection effects de- 
crease with radius and lead to underestimates of the mass- 
to-light ratio inside r2oo- Open squares (Figure 20) show 
the mass-to-light ratios of these shells. Indeed, there is 
a general trend for lower mass-to-light ratios in shells at 
larger radii, but the uncertainties are quite large. Note 
that the enclosed mass-to-light profiles M{< r)/L{< Rp) 
are weighted by mass and light and therefore differ from 
the (unweighted) profiles of mass-to- light ratios in individ- 
ual shells. For some clusters (e.g.. Coma and A119), the 
mass-to-light ratios of shells at large radii agree well with 
the enclosed mass-to- light profile, while for others (e.g., 
A496 and A576) the mass-to-light ratios in shells at large 
radii are significantly lower than the enclosed mass-to-light 
profile. This variety is likely due to projection effects, 
namely the large variety in the appearance of the caus- 
tic pattern for an individual cluster viewed from different 
lines of sight (D99). This variety suggests that the shapes 
of individual mass-to-light profiles should not be taken too 
seriously; however, the average mass-to-light profile should 
be unbiased (D99). 

X-ray mass estimates provide an independent check of 
our mass-to-light ratios within r^oo- As shown in Paper 
I, the caustic mass profiles evaluated at r^oo agree quite 
well with X-ray mass estimates from the mass-temperature 
relation. Thus, it is no surprise that the mass-to-light ra- 
tios calculated from the X-ray mass M500 and the lumi- 
nosity evaluated at L500 yield similar values. Note that, 
again, the luminosity estimates include all galaxies pro- 
jected within i?5oo- An NFW profile with concentration 
c=5 appropriate for clusters (Navarro et al. 1997) would 
have a deprojected mass-to-light ratio {M/L){< rsoo) ~ 
1.3M5oo/i(< -R500) (a higher concentration of c=20) re- 
duces this factor from 1.3 to 1.2). We show both the pro- 
jected and deprojected estimates as stars in Figure 20. 
Note that {M/L){< rsoo) « {M/L){< R500); these points 
may be compared with the profiles of mass in cylinders di- 
vided by light in cylinders. The mean deprojected mass-to- 
light ratio is {M/L){< rsoo) = 78±7(statistical)/iMQ/L0. 
This result agrees with the deprojected mass-to-light ra- 
tio {M/L){< r2oo) = 88 ± 9/1 (statistical) at r2oo taken 
from the caustic mass profiles assuming a correction of 
1.25 appropriate for a c=5 NFW halo. The agreement be- 
tween the mass-to-light ratios within rsoo and within r2oo 
siiggcsts that mass-to-light ratios are reasonably constant 
throughout the virial region of a cluster (see also §6.2). 
That is, the observed decrease in mass-to-light ratios with 
radius is not monotonic but may begin only at roughly 

Our results are in excellent agreement with Lin et al. 
(2003, hereafter L03), who find a mass-to-hght ratio 
M/Lk, = 76 ± 4/1 (statistical) at rsoo for hot [Tx > 
3.7keV) clusters using X-ray masses and 2MASS photom- 
etry from a larger cluster sample. We include only galax- 
ies within the caustics in our luminosity estimates; L03 
use statistical background subtraction to correct their lu- 



minosity estimates. The close agreement shows that the 
methods LOS use to subtract background galaxies do not 
introduce a bias in the luminosity estimates. 

Previously, we have used the caustic technique to calcu- 
late mass-to-light profiles in R band for A576 (Rines et al. 
2000) and in Kg band for Coma (Rines et al. 2001a). In 
A576, we found a steeply decreasing mass-to- light profile 
in R band. In Coma, we found a fiat profile but noted 
that the systematic effects allowed for a decreasing profile. 
The results we derive; lic^rc^ for these clusters are consistent 
with these earlier determinations. 

Other investigators have applied Jeans analysis to en- 
semble clusters to test for variations in the mass-to-light 
ratio. This effort is complicated by the fact that one needs 
to assume an orbital distribution to measure variations in 
M/L. Carlberg et al. (1997) and van der Marel et al. 
(2000) find that light traces mass in the CNOCl ensemble 
cluster (composed of massive clusters at z = 0.2 — 0.5) 
to a radius of 2r2oo- Biviano & Girardi (2003) construct 
an ensemble cluster from poor clusters in the 2dFGRS. 
They find that the ratio of the mass density to the galaxy 
number density decreases with radius to 2r2oo, similar to 
our result for the CAIRNS clusters. When they exclude 
late-type galaxies from the galaxy number density, the ra- 
tio is roughly constant. Katgert et al. (2004) construct 
an ensemble cluster from the ESO (European Southern 
Observatory) Nearby Abell Cluster Survey and find that 
the mass-to-light ratio decreases with radius in the range 
0.2 — 1.5r2oo, although the mass-to-light ratio is roughly 
constant when late-type galaxies are excluded. 

Weak lensing provides an independent estimate of mass- 
to-light ratios on large scales that does not depend on the 
dynamical state of the system. Kaiser et al. (2004) and 
Gray et al. (2002) estimate the mass-to-light ratios of su- 
perclusters with weak lensing. Kaiser et al. (2004) find 
M/Lb = 280 ± AOh for light in early-type galaxies. As- 
suming a typical early-type color of i? — Ks = 3.7 (Jarrett 
2000), this value corresponds to M/Lk^ w 64 ± 9h. In- 
cluding late-type galaxies would decrease this ratio. Wil- 
son et al. (2001) find similar results from weak lensing 
in blank fields. Gray et al. (2002) find M/Lb ^ 200h 
(early- type light only) for individual clusters; when they 
cross-correlate mass and light they find M/Lb = 130/i 
(early- type light only), but they caution that there are 
many systematic uncertainties in this estimate. Recently, 
Kneib et al. (2003) used weak lensing to estimate the mass 
profile of CL0024-hl654 to a radius of 3.25/i-^Mpc. Kneib 
et al. (2003) conclude that the K band mass-to-light ra- 
tio is roughly constant on these scales. Assuming passive 
evolution, their mass-to- light ratio corresponds to 65 ± 9h 
(74 ± lOh for red sequence galaxies only) at z = 0, inter- 
mediate between our estimates of the mass-to-light ratio 
inside and outside r2oo- 

Bahcall et al. (2000) use simulations to show that clus- 
ter mass-to-light ratios in B band exceed the global value 
due to the older, less luminous stellar populations found 
in clusters. Cluster mass-to-light ratios measured in Kg 
band should then be much closer to the global value be- 
cause Ks band light has a much weaker dependence on 
stellar population ages (e.g., Bell & de Jong 2001). If stel- 
lar populations are the primary cause of the decreasing 
mass-to-light profiles in the simulations of Bahcall et al. 
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(2000), the CAIRNS clusters should have roughly flat Ks 
band mass-to-light profiles. Thus, the similarity of the 
decreasing Kg band mass-to-light profiles of the CAIRNS 
clusters to the simulations of Bahcall et al. (2000) shows 
that the decreasing profiles in their simulations may not 
result primarily from differences in the stellar populations 
but from differences in the relative distribution of dark 
matter and galaxies. 

The CAIRNS sample is unique in both the completeness 
of the individual cluster catalogs and in the near-infrared 
digital photometry used to avoid stellar population effects. 
The mass-to-light ratios of the CAIRNS clusters decrease 
with radius and the mass-to-light ratios inside the virial 
regions agree with other estimates at optical and near- 
infrared wavelengths (see also Rines et al. 2000, 2001a). 
The decreasing mass-to-light profiles are consistent with 
results from other cluster studies. We discuss these re- 
sults in more detail in §7. 

6. PROPERTIES OF THE VIRIAL REGIONS 

6.1. The Halo Occupation Distribution 

The halo occupation distribution (see the review by 
Cooray & Shcth 2002) is an important input for convert- 
ing the results of numerical simulations into observablcs 
(e.g., Peacock k Smith 2000; Berlind & Weinberg 2002; 
Berlind et al. 2003, and references therein). The simplest 
prediction is that the number of galaxies formed is directly 
proportional to the available baryonic mass. Thus, the 
number of galaxies (brighter than some minimum mass 
or luminosity) contained in a halo of mass M is given by 
N (X M (i.e., the efficiency of galaxy formation is a uni- 
versal constant for sufficiently massive haloes). If galaxy 
formation is more efficient in the most massive haloes, then 
the relation might be A^ oc with /i > 1. Conversely, 
if galaxy formation is less efficient in massive haloes (e.g., 
if the gas is heated by the halo potential and is unable 
to collapse into galaxies) or galaxy disruption is more ef- 
ficient (e.g., dynamical friction and tidal stripping), then 
the relation might be A' oc M** with /z < 1. Models of 
the halo occupation distribution suggest that, for cluster 
mass halos, the relation is close to a power law with slope 
/X < 1. Semi-analytic models predict /i ^ 0.8 — 0.9 (Sheth 
& Diaferio 2001; Berlind et al. 2003). A smoothed particle 
hydrodynamics simulation of a ACDM cosmological model 
predicts halo occupation distributions with /i ^ 0.56 — 0.74 
for cluster mass halos, similar to the values for a different 
set of semi-analytic models (Berlind et al. 2003). Springel 
& Hernquist (2003) show that numerical simulations pre- 
dict suppression of galaxy formation in the most massive 
halos because gas cannot cool and collapse into galaxies. 

One of the few previous determinations of the relation 
between virial masses and galaxy numbers is that of Mari- 
noni & Hudson (2002), who compute masses and (blue) lu- 
minosities of virialized objects in the Nearby Optical Cata- 
log. Marinoni & Hudson (2002) find N cx M"-55±" "3, simi- 
lar to the semi-analytic models. Kochanek et al. (2003) use 
a constrained numerical simulation of 2MASS to develop a 
matched filter algorithm to study cluster properties in the 
2MASS catalog and heterogeneous auxiliary observations 
from the literature (e.g., redshifts and X-ray properties). 
Their best-fit relation between cluster mass and number 
of members is A' oc Mi-i°±°-°9. Pisani et al. (2003) find 



A^ cx 7i^o 'i'o±o.o4 jj-^ sample of groups, although mass es- 
timates of groups are very uncertain. Recently, Lin et al. 
(2004) analyzed the halo occupation distribution for clus- 
ters with 2MASS photometry and X-ray mass estimates. 
They find A^ cx Af" **'*±" °'', steeper than Marinoni & Hud- 
son (2002) but still in reasonable agreement with models. 

We can constrain the halo occupation distribution with 
the CAIRNS clusters, which cover roughly an order of 
magnitude in mass and have both accurate photometry 
and complete spectroscopy. Our resiilts have the advan- 
tages of uniform sky coverage, greater redshift complete- 
ness, and galaxy selection at near-infrared wavelengths, 
which is a better tracer of stellar mass and suffers less 
dust extinction than blue light. Conveniently, the mag- 
nitude limit we adopt {Mk, < M^r -\- 1) is very similar 
to the luminosity threshold used in Berlind & Weinberg 

(2002) and one of the thresholds used in Berlind et al. 

(2003) , L > 0.5L*. Figure 21 shows the number of galax- 
ies A'200 projected within i?200 versus M200, the mass of 
the halo. We do not attempt to deproject the number 
density profiles to obtain a deprojected estimate of A^200 
because the number density profiles are too noisy. If all 
haloes have similar concentrations, then the fraction of 
interlopers should be constant with mass. If the halo 
concentration decreases with mass (as expected for NFW 
models), then the fraction of interlopers should increase 
with mass. In this case, the fit to fi would be an over- 
estimate. The bisector of the two ordinary least-squares 
fits (Feigelson & Babu 1992) yields A^aoo oc M^o^°=^°•°^ 
3.3(7 shallower than a linear relation (shown by a dashed 
line in Figure 21). This result is not driven by A 194, 
the least massive cluster; excluding this cluster yields a 
least squares fit A-200 oc M°5^^±°-i^. This result agrees 
well with previous determinations as well as with expecta- 
tions from semi-analytic models for galaxy formation (e.g., 
Kauffmann et al. 1999a; Sheth & Diaferio 2001; Marinoni 
& Hudson 2002; Berlind et al. 2003; Pisani et al. 2003; Lin 
et al. 2004). We speculate that the significant difference 
from Kochanek et al. (2003) is due to the systematic un- 
certainties from the process of matching their simulation 
to the observations. Kochanek et al. (2003) use a matched 
filter algorithm which is finely tuned to reproduce the ex- 
pected properties of clusters based on simulations (where 
galaxies trace the dark matter distribution). Systematic 
effects can arise both from mismatches in the assumed and 
true cosmology and recipes for galaxy formation as well as 
unknown systematics in the heterogeneous auxiliary ob- 
servations. 

The comparison with Lin et al. (2004) is especially in- 
teresting because both datasets use 2MASS photometry. 
Lin et al. (2004) use a much larger sample of clusters but 
they use only statistical background subtraction whereas 
we study fewer clusters but use complete spectroscopic in- 
formation to assign cluster membership. A detailed com- 
parison of these two methods would be instructive but it 
lies beyond the scope of this paper. In particular, there 
are few clusters in both samples, so cluster-to-cluster vari- 
ations could significantly affect the comparisons. We refer 
the reader to Lin et al. (2004) for an excellent discussion 
of the physical significance of a non-linear HOF as well as 
the observational implications for clusters. 
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6.2. Mass Dependence of Mass-to-Light Ratios 

Figure 22 shows the mass-to-Ught ratio within r2oo ver- 
sus M200 for the CAIRNS clusters. The scatter is large, 
but the CAIRNS clusters show an increase in M/L (evalu- 
ated at r2oo) with increasing mass. Lin et al. (2003, here- 
after L03) found a similar correlation between X-ray mass 
and near- infrared mass-to- light ratios; more massive clus- 
ters have larger mass-to-light ratios with a best-fit relation 

-^500 ^0.31±0.09 



(M/LkJ(< rsoo) = (67±4)/i( 



2.1 X 10"/i-iAffr 



(7) 

Figure 22 shows this relation assuming that the mass-to- 
light ratios within r5oo and r2oo are identical, M200 ~ 
3M500, and multiplying by 0.8 to convert to the spheres- 
by-cylinders measured here. The CAIRNS clusters follow 
this relation quite closely, showing that the mass estimator 
used (X-ray versus virial/caustic mass) does not affect the 
correlation of mass-to-light ratio with cluster mass. The 
close agreement also demonstrates that cluster mass-to- 
light ratios do not change dramatically between rsoo and 
^■200 (see also § 5). 

A compilation of virial masses and luminosities by Gi- 
rardi et al. (2000) yields (M/Lb,) oc mO-I^-O-^^ Simi- 
larly, Bahcall & Comcrford (2002) use a heterogeneous 
catalog to derive a dependence of (optical) M/L on X- 
ray temperature which they attribute to differences in 
the ages of stellar populations. One can convert their 
relation into a M/L — M relation with the X-ray mass- 
temperature relation (Finoguenov et al. 2001). Specif- 
ically, (M/L) oc t0-30±0 08 and M oc Ti-64±0.04 ^-^^^ 

(M/L)|r2oo oc M2oo^^°°^, sUghtly shallower than but in 
agreement with the LOS relation and the CAIRNS clusters 
(Figure 22). Note, however, that the CAIRNS relation has 
little dependence on the ages of the stellar populations, 
counter to the conclusion of Bahcall & Comerford (2002). 
If differences in stellar populations produce the M/ L — M 
relation, the slope of the relation should be steeper for 
M/Lb than for M/Lk,- 

In contrast, Kochanek et al. (2003) find that mass-to- 
light ratios are smaller in more massive clusters; they 
find a best-fit relation of {M/LkJ{< ?'20o) = 116 ± 
46(M2oo/10^^/i~^Mq)-o io±o o9. We multiply this rela- 
tion by 0.8 to convert to the spheres-by-cylinders mea- 
sured here (Figure 22). The CAIRNS clusters fofiow the 
relation foimd by LOS, Girardi et al. (2000), and Bahcall 
et al. (2000), and exclude the relation of Kochanek et al. 
(2003). The disagreement with Kochanek et al. (2003) 
is perhaps not surprising given the disagreement between 
their halo occupation function and that of the CAIRNS 
clusters found in the previous section. 

Figure 2S shows the relation between L200J and N200, 
the number of bright galaxies {Mk, < —22.77 -|- 51og/i) 
projected within i?200- The bisector of the ordinary least- 
squares fits is A^200 oc L^.Q^^^'''''^'^ , consistent with a slope 
of unity. This result underscores the result of §S.l; the 
cluster-to-cluster variations in the LF are small. Galaxy 
formation is suppressed (and/or that the efficiency of 
galaxy disruption is enhanced) in massive clusters, with 
greater suppression for more massive clusters. The corre- 
lation of mass-to-light ratio with mass is then a natural 
byproduct of the correlation of M2Q0/N2QQ with mass and 
a universal LF. 



These results are consistent with the decreasing mass- 
to-light profiles found in §5. These decreasing profiles 
imply that cluster infall regions, which contain galaxies 
formed in environments with smaller virial temperatures 
than galaxies in the virial regions, have smaller mass-to- 
light ratios. The presence of X-ray groups in cluster infall 
regions (Rines et al. 2001b, 2002) demonstrates the over- 
lap between cluster infall regions and low-mass clusters. 
In §6.1, we show that the number of galaxies within -R200 
increases more slowly than the cluster mass. Thcisc results 
all suggest that the efficiency of galaxy formation is sup- 
pressed (see the numerical simulations by Springel & Hern- 
quist 2003) and/or that the efficiency of galaxy disruption 
is enhanced (see the numerical simulations by Kravtsov 
& Klypin 1999; Colin et al. 1999) in environments with 
larger virial temperatures. In the latter case, the contri- 
bution of intracluster stars to the cluster light budget can 
be substantial (5-50%, see §7.4). Neglecting this contri- 
bution (the normal procedure and the one adopted here) 
may result in a severe underestimate of the total light in 
the cluster. 



7. DISCUSSION 

7.1. Predicting Mass Profiles From the Galaxy 
Distributions 

Because the caustic technique is relatively new, we test 
the consistency of our results with velocity dispersion pro- 
files (VDPs), a more traditional tool of galactic dynam- 
ics. Here we predict the mass profiles based on the ob- 
served distribution of galaxies assuming that they trace 
the mass. If the radial variations in the mass-to-light ra- 
tio (§5) are real, then the mass profiles calculated from the 
galaxy distributions (either number density or luminosity 
density) assuming a constant mass-to-light ratio should 
differ demonstrably from those in Paper I. We reproduce 
the VDPs from Paper I in Figure 24 (dash-dotted lines). 

In §4.1 above, we fit the surface number density profiles 
of the bright galaxy distribution. The scale and normaliza- 
tion of the profiles predict the shape of the velocity disper- 
sion profiles (VDPs) of the clusters under the assumption 
of isotropic; orbits and a (globally) constant ratio M /N of 
mass to number of bright galaxies (M^^ < — 22.77-|-51og/i). 
We calculate the value of this ratio within r2oo and find 
iM/N)\r2Qa = (7.8 ±0.8) x lO^^/j-^Megalaxy-^ Figure 
24 displays these predicted VDPs (solid lines) along with 
the observed VDPs and the VDPs of the best- fit Hernquist 
mass profiles from Paper I (dash-dotted lines). As noted in 
Paper I, the VDPs predicted by the caustic mass profiles 
agree well with the observed VDPs. The VDPs predicted 
from the surface density profiles, however, do not agree 
with the observed VDPs, especially in Coma and A576, 
where the predicted VDPs lie substantially below the ob- 
servations. 

This disagreement may result from anisotropic orbits or 
from cluster-to-cluster variations in M/N (we find such 
variations in §6.1). Anisotropic orbits would affect both 
the shapes and normalizations of the predicted VDPs, but 
variations in M/N only affect the normalizations (assum- 
ing that M/N is constant with radius in each cluster). Fig- 
ure 23 shows that {M/N)\r2oo increases with M2oo- Cor- 
recting for this trend reconciles some of the differences in 
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Figure 24, i.e., the most discrepant clusters are those with 
the highest masses. Substituting a higher {M/N)\r2oo in- 
creases the amphtude of the predicted VDPs and brings 
the predicted and observed VDPs into better agreement. 
Outside r2oo, galaxies are not relaxed. At these radii, 
VDPs do not necessarily contain information about the 
orbital distribution. Thus, the VDPs outside r2oo should 
not be considered strong constraints on the mass and or- 
bital distributions. 

The most straightforward predictions of VDPs based on 
the assumption that light traces mass disagree with the 
observed VDPs. In contrast, the VDPs predicted by the 
caustic mass profiles agree well with the observed VDPs. 
Thus, the decrease in the efficiency of galaxy formation 
(and/or the increase in the efficiency of galaxy disruption) 
for haloes with larger virial temperatures is not an arti- 
fact of the caustic technique. Velocity dispersion profiles, 
a more traditional tool of galactic dynamics, also indicate a 
discrepancy between the distribution of galaxies and mass 
(although subject to possible biases from the orbital dis- 
tribution and/or the dynamical state of galaxies outside 
?'2oo)- It is interesting that Katgert et al. (2004) also find 
a decreasing mass-to-light profile for an ensemble cluster 
using Jeans analysis to compute the mass profile. Their 
results strengthen our conclusion that the decreasing mass- 
to-light profiles are physical effects. 

7.2. Morphological Gradients 

Because of the well-known correlation between morphol- 
ogy and density, we expect larger fractions of late-type 
galaxies with increasing clustrocentric radius (decreasing 
density). If the LFs of early- type and late- type galaxies 
differ significantly, the total LF should vary with clustro- 
centric radius. That is, the LF in cluster centers should 
closely resemble the early-type LF, whereas at larger radii 
it should resemble the late-type LF. KOI separate the LF 
into early-type and late-type LFs and find that fits to 
Schechter functions yield a brighter M|-^ for the early-type 
LF; the faint-end slope is slightly shallower for the late- 
type LF (but see Bell et al. 2003, who find that 2MASS 
misses many blue low surface brightness galaxies present 
in SDSS). 

Because our spectroscopic surveys extend to fixed abso- 
lute magnitudes, the correction for light in galaxies fainter 
than our limiting magnitudes changes with radius. Adopt- 
ing the type-dependent LFs of KOI, the correction be- 
comes larger with radius provided the limiting magni- 
tude is Mk^ ^ (note that the type-dependent LFs 
are only constrained for MK^,iso < —20.5). That is, a 
magnitude-limited survey misses more light at large clus- 
trocentric radii. If the fraction of early- type galaxies 
changes from 1 to (a huge overestimate), the correc- 
tion changes by ^20% (~10%) for a magnitude limit of 
-^ifs,iso=-22.5 (-21.5), approximately the magnitude lim- 
its for the CAIRNS clusters. Correcting for this effect in 
§5 would add even more light to the cluster outskirts and 
lead to more steeply decreasing profiles. Thus, the de- 
creasing mass-to-light profiles in §5 cannot be explained 
by the morphology- density relation. 

7.3. Stellar Populations and Correcting for Faint 
Galaxies 



Profiles of the ratio of dark matter to stellar mass can 
be used both to estimate fim and to constrain prescrip- 
tions for galaxy formation. If iiT^-band light traces stellar 
mass exactly, the results of §5 indicate that the efficiency 
of star formation is reduced in dense cluster environments. 
However, the properties of galaxies change rapidly with in- 
creasing distance from cluster centers (Balogh et al. 2004, 
and references therein). In particular, the stellar mass- to- 
light ratio is smaller in late-type galaxies than in early-type 
galaxies by up to a factor of two. Late-type galaxies are 
much more common in the field than in clusters. Thus, 
the mean stellar mass-to-light ratio should decrease with 
radius. 

Because stellar populations are younger at larger clus- 
trocentric radii, mass-to-optical-light profiles might de- 
crease with radius (Bahcall et al. 2000) even if the ratio of 
gravitational mass to stellar mass is constant. Thus, the 
total-to-stcllar mass profiles of the CAIRNS clusters may 
decrease less steeply than the mass-to-light profiles. Here, 
we test for radial gradients in the stellar mass-to-light ra- 
tio in A576 and estimate the magnitude of this effect in 
Kg band. 

7.3.1. A Test in A576 and the Importance of Faint 

Galaxies 

We showed above (§3.3) that there are no obvious gra- 
dients in the R — Kg colors of galaxies in A576. We test 
for gradients in the stellar mass-to-light ratio directly by 
comparing the mass-to- light profiles in an optical band (R) 
and near-infrared bands. If there were a significant gradi- 
ent in stellar mass-to-light ratios, the near-infrared profile 
would be flatter than at optical wavelengths. 

Indeed, the Kg band mass-to-light profile (thick solid 
line in Figure 25) decreases more slowly than the R band 
profile (dash-dotted line). The cumulative mass- to-light 
ratio decreases by a factor of ~2 in R band and by a fac- 
tor of ^1.4 in Kg band. This result suggests that the 
effect of star formation gradients on mass-to-light profiles 
is significant in the R band as well as in the B band (e.g., 
Bahcall et al. 2000). However, there is no obvious change 
in the average R — Kg color with radius in A576 (§3.3). 
Thus, the steeper decrease in the cumulative mass-to- light 
profile in the R band is not readily explained by a simple 
color gradient. 

Another explanation for the difference in R and Kg band 
M / L profiles is the corrections for faint galaxies without 
redshifts. The R band catalog has complete spectroscopy 
to R=16.5 and complete photometry to R=18.0. Rines 
et al. (2000) used several techniques to estimate the (as- 
sumed constant) flux surface density contributed by back- 
ground galaxies. Here, by contrast, we correct for faint 
galaxies by assuming a universal luminosity function in all 
environments. Under this assumption, the total luminos- 
ity in all galaxies is simply a constant factor multiplied by 
the luminosity contained in bright galaxies. 

Because 2MASS is a shallow survey, it is difficult to es- 
timate magnitudes (and hence number counts) of galaxies 
fainter than our spectroscopic completeness limit (these 
number counts are necessary to make a background cor- 
rection similar to Rines et al. 2000). However, it is interest- 
ing that the 2MASS galaxy counts indicate a steeply rising 
LF in the outskirts of A576 fainter than the spectroscopic 
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completeness limit (Figure 4). The number counts there- 
fore suggest that applying a constant background subtrac- 
tion to the 2MASS data would lead to better agreement 
between the two mass-to-light profiles. The presence of 
a background group or cluster behind A576 (Rines et al. 
2000) shows that the true background is non-uniform. 

Fortunately, it is straightforward to apply the assump- 
tion of a universal luminosity function to the mass-to- 
light profile in the R band. Under this assumption, the 
two mass-to-light profiles arc in quantitative agreement; 
the decrease in the mass-to-light ratio between the inner 
1 h~^Mpc and the outskirts is a factor of ~1.4 at both 
wavelengths. Figure 25 shows that the shapes of the mass- 
to-light profiles (thick solid and dashed lines) agree well 
under this assumption (except at Rp < 0.5 h~^M.pc where 
the photometry of two bright galaxies is uncertain; see 
§2.2). Thus, the apparent disagreement between the R- 
band and fCs-band shapes of the cumulative mass-to-light 
profiles in A576 is not due to radial gradients in the stellar 
mass-to-light ratio but is simply a result of using different 
methods to account for the luminosity contributed by faint 
galaxies. Using consistent methods produces both qualita- 
tive and quantitative agreement in the mass-to-light pro- 
files calculated at optical and near-infrared wavelengths. 
Because the same caustic mass profile is used for both 
wavelengths, this comparison tests the relative shapes of 
the light profiles calculated at different wavelengths. The 
agreement between the two profiles should therefore gen- 
eralize to all clusters with similar galaxy populations; the 
fact that A576 has one of the most strongly peaked mass- 
to-light profiles (Figure 22) should not affect this gener- 
alization. This result demonstrates the importance of us- 
ing consistent corrections for faint galaxies when compar- 
ing mass-to-light ratios or profiles at different wavelengths 
and/or for different clusters. This result also suggests that 
a complete census of cluster light requires deep, complete 
spectroscopy. 

7.3.2. Estimating the Total-to- Stellar Mass Profiles 

Although we showed above (§3) that there is no obvious 
evidence from near-infrared photometry or i? — Ks colors 
for dramatic changes in the stellar populations with clus- 
trocentric radius, the degeneracy between age and metal- 
licity and the weak dependence of near-infrared colors on 
these properties might obscure a real gradient. Here we 
estimate the potential size of this effect. 

Because the galaxies we sample are relatively bright, 
the color-magnitude relation implies that these galaxies 
should have red colors. We estimate the fraction of light in 
early-type galaxies at large radii using the type-dependent 
LFs of KOI. Using their estimates of the type-dependent 
LFs, early-type galaxies contribute roughly half of the to- 
tal hght in bright galaxies {Mr, < -22.77 -h 51og/i) av- 
eraged over all environments. The Kg band mass-to- light 
ratio is a factor of 1.8 ± 0.3 times larger in virial regions 
than in infall regions (§5). Thus, the ratio of total matter 
to stellar matter (in galaxies) could be roughly constant on 
scales up to ^ 10 ft'-'^Mpc if the Kg band stellar mass-to- 
light ratio in early-type galaxies is ^2.6 times larger than 
in late-type galaxies. Bell & de Jong (2001) and Bell et al. 
(2003) find that the stellar mass-to-light ratio measured in 
Ks band varies by no more than a factor of 2 over a wide 



range of star formation histories. Such a large difference 
would probably produce significant R—Kg color gradients, 
contradicting Figure 13. The age-metallicity degeneracy 
and/or complicated star formation histories could conceiv- 
ably mask these gradients, but these effects are generally 
small at Kg band. 

A second method of quantifying the changes in stellar 
mass-to-light ratios is to use the relation between stellar 
mass-to- light ratio and galaxy color (Bell & de Jong 2001; 
Bell et al. 2003). The range of galaxy colors in the SDSS is 
0-4 ^ g — r < 1.0, although little stellar mass is contained 
in the bluest galaxies (e.g., Kauffmann et al. 2003). If the 
average galaxy g — r color changed from 0.9 in cluster cen- 
ters to 0.4 in the outskirts (an extreme assumption), the 
corresponding change in stellar mass-to-light ratio is only 
a factor of 1.3. A more realistic estimate is that early- 
type red galaxies comprise roughly half the light in bright 
galaxies. Then, the average Kg band stellar mass-to-light 
ratio in cluster outskirts is at most 20% smaller than in 
cluster centers. Thus, gradients in stellar populations do 
not account for the radially decreasing mass-to-light ratios 
of the CAIRNS clusters. 

7.4. Intracluster Light 

Another explanation for the decreasing mass-to-light 
profiles is that we do not account for light outside of galax- 
ies. The existence of intracluster red giant branch stars 
(Durrell et al. 2002), planetary nebulae (Ciardullo et al. 
1998; Feldmeier et al. 1998; Durrell et al. 2002; Arnaboldi 
ct al. 2003; Feldmeier et al. 2003), globular clusters (West 
et al. 1995; Jordan et al. 2003), diffuse light (Zwicky 1951; 
Melnick et al. 1977; Uson et al. 1991; Bernstein et al. 1995; 
Gregg & West 1998; Trentham & Mobasher 1998; Gonzalez 
ct al. 2000; Feldmeier et al. 2002). and supernovae (Smith 
1981; Gal- Yam et al. 2003) not associated with individual 
galaxies all suggest that stars are stripped from cluster 
galaxies and form diffuse intracluster light (Moore et al. 
1999; Gnedin 2003). Numerical simulations of clusters 
in ACDM cosmologies show that processes such as tidal 
stripping and dynamical friction disrupt cluster galaxies 
(Kravtsov & Klypin 1999; Colin et al. 1999). 

Surveys of the above tracers of intracluster stars indicate 
that intracluster light constitutes ^5-50% of the total light 
in the virial regions. The decreasing mass-to-light profiles 
found here may be flat if intracluster light is taken into 
account. Thus, the most secure conclusion we can draw 
is that the number of galaxies per unit mass is smaller 
in cluster virial regions than in infall regions. The star 
formation efficiency could be constant in all environments 
with the observed dependence resulting from more efficient 
galaxy disruption in environments with larger virial tem- 
peratures. A prediction of this scenario is that the fraction 
of intracluster light should increase with cluster mass. A 
similar trend has recently been noted in simulated clusters 
(Murante et al. 2004). 

7.5. Uncertainties in the Caustic Mass Profiles 

D99 used numerical simulations to investigate the sys- 
tematic uncertainties in the caustic technique. D99 found 
that the uncertainties in individual cluster profiles are 
large (~ 50%) but unbiased. We extend this work to ob- 
servations in Paper I. Surprisingly, the contrast between 
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the caustic envelope and the background is larger in the 
CAIRNS clusters than in the simulations of D99. This 
difference may indicate a mismatch between the cosmo- 
logical model used in D99 (standard ACDM) and the true 
model and/or deficiencies in the recipe for star formation 
and galaxy formation used in the simulations. 

D99 anaylze only Coma-size clusters. It is possible that 
the caustic technique is less accurate and/or biased for less 
massive halos. Further high-resolution simulations with 
different cosmological models and/or different recipes for 
star and galaxy formation may clarify this issue. 

The masses obtained with the caustic technique agree 
very well with virial masses and X-ray estimates at small 
radii (Paper I). Thus, the mass-to-light ratios at these radii 
are reasonably secure (if corrected for projection effects). 
The caustic mass profiles in Paper I agree very well with 
both NFW and Hernquist models, with each model pro- 
viding a slightly better fit on roughly half of the clusters. 
The NFW profile predicts more mass at large radii than 
the Hernquist profile, and it produces better fits to halos 
in CDM simulations (Navarro et al. 1997). The extrapola- 
tion of the NFW profile beyond the virial radius provides 
a reasonable description of clusters in simulations (Tasit- 
siomi et al. 2004). We therefore calculate the mass-to- light 
ratios assuming that the NFW profiles inside r2oo extend 
to Tt and that the caustic diagrams indicate cluster /infall 
region membership. This calculation yields results similar 
to those foimd using the caustic mass profiles. 

Caustics are a good but not perfect indicator of clus- 
ter/infall region membership. Galaxies outside the caus- 
tics are outside the infall region, but there may be inter- 
lopers in the caustic diagram. In cluster cores, only ~1% 
of galaxies arc interlopers (van der Marel et al. 2000). The 
number of interlopers should increase at roughly the same 
rate as the area sampled. Indeed, detailed numerical sim- 
ulations of clusters indicate that the fraction of interlopers 
increases with radius (Cen 1997). To estimate this effect 
in Coma, we estimate the contribution of interlopers to 
be the background luminosity density times the volume 
within the caustics, which scales roughly with the area on 
the sky. The luminosity contained in possible interlop- 
ers in each radial bin is less than 10% of the luminosity 
in that bin. Thus, although the luminosity of interlopers 
could lead to an overestimate of the luminosity of infall 
region members at large radii, the overestimate is likely 
<10%. 

Finally, §7.1 shows that the observed VDPs (a better es- 
tablished tool of galactic dynamics) differ from those pre- 
dicted by the galaxy distributions under the assumption of 
isotropic orbits. The above discussion indicates that the 
decreasing mass-to-light profiles are probably not caused 
by (currently unknown) systematic effects in the caustic 
technique. 

7.6. Cosmological Implications 

Two groups have used the Second Incremental Data 
Release (2IDR) of 2MASS to measure the near-infrared 
galaxy luminosity fimction (K01,C01). Both groups find 
acceptable fits with the functional form proposed by 
Schechter (1976). Further, their estimates of the near- 
infrared luminosity density and the best-fit parameters of 
the LF agree within the uncertainties. Assuming that 



the cluster luminosity function is identical to the field 
luminosity function of COl (calculated using 2MASS Kg 
extrapolated magnitudes, the same magnitude definition 
used here), the average mass-to- light ratio within the 
turnaroimd radius {M/ LK,)tot implies flm = 0.18 ±0.04 
(statistical) . If the global value of M/Lk, is closer to the 
value in cluster infall regions than the value in cluster virial 
regions, the best estimate of Q.m is from {M/LKjinf (the 
average mass-to- light ratio between r2oo and rt), which 
yields Qm ~ 0.13 ± 0.03 (statistical). Note that these es- 
timates become ~20% smaller if we adopt the COl 2IDR 
Kron magnitude LF with a -0.20 magnitude adjustment 
for converting 2IDR Kron magnitudes to 2MASS extrapo- 
lated magnitudes (the best-fit LF parameters for the Kron 
LF have a fainter and a shallower faint-end slope). 

We explicitly use the luminosity function of COl (for 
extrapolated magnitudes) to estimate the completeness 
corrections for faint galaxies. This constraint means that 
our results are independent of the faint-end slope of the 
luminosity function in clusters. Provided the luminosity 
functions are similar at the bright end (as shown in §3.1), 
the estimate of 17^ is independent of the properties of 
dwarf galaxies. Similarly, any unusual systematic effects 
in the measured photometric properties of the galaxies 
are present in both our sample and that of COl. Thus, 
any such effects should cancel out in the estimate of 
In particular, 2MASS misses faint, low surface brightness 
(LSB) galaxies (Andreon 2002a; Bell et al. 2003). Be- 
cause these galaxies are missing in both the CAIRNS cat- 
alogs and in the estimates of the field LF, the omission of 
these galaxies leads to an overestimate of M / Lk^ but does 
not affect the estimate of ilm- If these LSB galaxies were 
substantially more numerous in low-density environments 
than in cluster environments, a bias could result, but An- 
dreon (2002a) shows that these LSB galaxies are present 
in clusters. Furthermore, most of these LSB galaxies are 
fainter than the portion of the LF sampled in the CAIRNS 
2MASS catalogs {Mk, < -22). 

Some investigators suggest that the local universe is sub- 
stantially underdense with respect to the global average 
density (e.g., Busswcll et al. 2003; Frith et al. 2003, and 
references therein). In particular, galaxy number counts 
indicate that the region surveyed by the 2dFGRS is signifi- 
cantly underdense (Frith et al. 2003) . The presence of such 
an underdensity obviously has important implications for 
estimating ilm using the mass-to-light ratio. 

Wright (2001) suggests that the near-infrared luminos- 
ity density estimated by COl is a factor of 2.3 smaller 
than the value obtained by extrapolating the z band Sloan 
Digital Sky Survey LF using typical spiral galaxy colors. 
However, Blanton et al. (2003) recently released a cor- 
rected version of the SDSS LF (including evolutionary K- 
corrections) which yields different LF parameters and a z 
band luminosity density smaller by a factor of 1.29. From 
a comparison of galaxies in both SDSS and 2MASS, Blan- 
ton et al. (2003) find that the mean difference between 
the °-^i SDSS band (the notation means that the band- 
pass is the rest- frame bandpass of a galaxy at z = 0.1 as 
observed in the SDSS i band, i.e., "'^z is slightly blueward 
of ° the observed bandpass) and the 2MASS Kg band 
is^'^i — i<'s~2.52. The ^'^i band luminosity density can 
then be extrapolated to Kg band (using Mqo.h = 4.58 and 
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Mq^k, = 3.39) to obtain j{Ks) ~ 7.22 x lO^/iioMpc"^. 
We thus obtain fi^ ~ 0.14 ± 0.05 using the total (virial 
plus infall region) mass-to- light ratio {M/LKjtot, and 
w 0.10±0.03 using the mass-to- light ratio (M/ Lk^ )inf 
in the infall region only. 

Huang et al. (2003), using a smaller but deeper sur- 
vey, suggests that the infrared luminosity density is signif- 
icantly larger (but see Bell et al. 2003, who note that this 
estimate ignores evolutionary corrections). Using his lumi- 
nosity density yields flm = 0.23 ±0.04 from M/LkJ< rt) 
and = 0.16 ± 0.03 from M/LK,{r200 ft)- Including 
evolution corrections reduces the estimate of flm by ~20% 
(Bell et al. 2003), yielding estimates similar to those for 
SDSS and 2dFGRS. Clearly, the normalization of the in- 
frared luminosity density is a significant source of uncer- 
tainty in using the cluster mass-to-light ratio to determine 
the matter density. 

Table 9 list several recent estimates of il„i from a variety 
of techniques. The estimates in Table 9 typically assume 
a flat universe dominated by dark energy. A detailed dis- 
cussion of the systematic uncertainties and potential biases 
in the various techniques lies outside the scope of this pa- 
per. In general, estimates of Qm from cluster abundances 
and dynamics and weak lensing yield low values; estimates 
from supernovae and from the combination of microwave 
background with large-scale structure yield higher values 
of flm- Our estimates are smaller than the currently pop- 
ular value of il„i ~ 0.27, but within the range of estimates 
from other techniques. It is curious that our estimates 
agree with other estimates based on mass-to-light ratios 
both inside and outside of clusters. 

However, we find a significantly smaller value than es- 
timates based on the cluster baryon fraction. In particu- 
lar, L03 (see also Mohr et al. 1999) calculate the baryon 
fraction within rsoo and estimate flm = 0.28 ± 0.03 (sta- 
tistical). Comparing this estimate with the mass-to-light 
ratios of hot clusters in their sample, they conclude that 
the mass-to- light ratio in hot clusters {kTx > 3.7keV) is 
a factor of 0.68 ± 0.10 smaller than the global value. This 
conclusion disagrees with our result that the mass-to-light 
inside r2oo is a factor of 1.8 ± 0.3 larger than the mass-to- 
light ratio outside r2oo (which should better approximate 
the global value). Note that departures from hydrostatic 
equilibrium of intracluster gas due to nonthermal pres- 
sure would aggravate this problem by decreasing the true 
cluster baryon fraction (e.g., Sadat & Blanchard 2001). 
At least two explanations may account for the discrep- 
ancy between L03 and the decreasing mass-to-light pro- 
files. First, if baryons in hot gas avoid cluster centers due 
to, e.g., shock heating, the baryon fraction within rsoo 
may be smaller than the global value (see the compari- 
son of many simulations by Frenk et al. 1999). Detailed 
observations with ROSAT and ASCA showed that the gas 
mass fraction increases with radius in some nearby clusters 
(David et al. 1995; Markevitch & VikhUnin 1997; Ettori & 
Fabian 1999), but little data exist beyond rsoo. Even with 
Chandra and XMM-Newton, warm/hot gas is presently not 
observable at large radii because its temperature and den- 
sity are too low. If the baryon fraction continues to in- 
crease outside rsoo, the baryon fraction within rsoo leads 
to an overestimate of flm- Thus, it may be possibile to 
reconcile these results, but at present, the baryon frac- 



tion outside rsoo remains unconstrained by observations. 
Second, cluster infall regions and low-mass clusters (which 
L03 and §6.2 show have smaller mass-to-light ratios than 
more massive clusters) may provide a uniquely favorable 
environment for star formation where the baryon density 
is high enough to encourage gravitational collapse but not 
so high that virial temperatures prevent collapse. Under 
this scenario, the mass-to-light profile of a cluster would 
peak in the center, decrease in the infall region, then rise 
again to the global value. 

The latter explanation is intriguing, and might lead to 
consistency with Turner (2002), who notes that flm = 0.33 
is signicantly larger than previous determinations based on 
the mass-to-light ratios in clusters and concludes that the 
difference results from variations in the mass-to-light ra- 
tio with environment. The mass-to- light profiles presented 
here disagree with this conclusion both qualitatively and 
quantitatively to densities as small as w Spc- The mass-to- 
light ratio decreases with radius, and there are no obvious 
systematic effects in the caustic technique that can rec- 
oncile these results. Thus, the scale dependence of the 
mass-to-light ratio on scales < 10 h~^Mpc cannot be the 
resolution of this profound problem. 

Some recent simulations suggest that galaxies form pref- 
erentially in overdense regions of the universe (Blanton 
et al. 1999; Ostriker et al. 2003). These simulations im- 
ply that the estimates of flm from the mass-to-light ratio 
in c;lustcr virial regions may underestimate the true value 
by a factor of ~1.25. We find instead that the mass-to- 
light ratio decreases in cluster infall regions. It is possible, 
however, that the global mass-to-light ratio is significantly 
higher than in cluster infall regions and low-mass clusters. 
Although somewhat arbitrary, such a scenario is consis- 
tent with all the constraints found in this paper and other 
investigations. Future studies of bulk flows are the most 
likely candidate to test this scenario. 

8. CONCLUSIONS 

We discuss some of the first estimates of radial variations 
in mass-to-light ratios on scales of l-10ft-~^Mpc using near- 
infrared photometry from 2MASS and mass profiles from 
the kinematics of inf ailing galaxies. Because cluster infall 
regions contain the transition from cluster galaxies to field 
galaxies (EUingson et al. 2001; Lewis et al. 2002; Gomez 
et al. 2003; Treu et al. 2003; Balogh et al. 2004, and refer- 
ences therein), mass-to- light ratios in infall regions should 
closely resemble the global value. 

To summarize our results: 

• Infall regions contain more bright galaxies (to a 
fixed absolute magnitude limit) than cluster virial 

regions. 

• The near-infrared luminosity functions for bright 
galaxies (Af^,, < -22 + 51og/i) in the CAIRNS 
cluster virial regions and infall regions do not 
differ significantly from the field galaxy luminosity 
function. Clusters contain an excess of extremely 
bright galaxies above the predictions of a Schechter 
function. 

• Optical-near-infrared colors in A576 show no 
radial dependence. This lack of a color gradient 
shows that the stellar populations do not change 
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dramatically with radius. It is likely that the mild 
gradients found in the photometric study of Goto 
et al. (2004) may be enhanced in optically selected 
samples as compared to near-infrared selected 
samples such as CAIRNS. 

Galaxies in cluster virial regions and infall regions 
exhibit a near-infrared color-magnitude relation 
with a shallower slope than at optical wavelengths. 

These galaxies also exhibit little scatter in J — Kg 
colors, indicating that the stellar populations 
are fairly homogeneous and that internal dust 
extinction and/or emission is important for only a 
few galaxies. 

Both the surface number density profiles and 
surface luminosity density profiles of CAIRNS 
members indicate that galaxies and stellar light 
are more extended than mass. 

Near-infrared mass-to-light ratios generally 
decrease with radius by a factor of 1.8 ± 0.3 in the 
infall regions of the CAIRNS clusters. This result 
agrees with previous results based on individual 
clusters and optical photometry. The presence of 
decreasing mass-to-light profiles even at Kg band 
suggests that the decrease is not due to changes in 
stellar populations. 

Near-infrared mass-to-light ratios calculated at 
r2oo using caustic mass estimates agree quite well 
with mass-to-light ratios calculated at rsoo from 
X-ray mass estimates. This agreement suggests 
that the decreasing mass-to-light profiles are not 
monotonic; the mass-to-light ratio is roughly 
constant inside r2oo- 

We derive some of the first constraints on the 
halo occupation function using cluster masses 
and near-infrared selected galaxy samples. The 
number of bright galaxies N2oa projected within 
i?200 increases as 

N200 oc M°o^°^° °^ significantly 
shallower than N200 oc M2oo- Earlier studies of 
the halo occupation distribution suggest that 
a halo occupation distribution shallower than 
A^200 oc M200 is necessary to reproduce the 
observed clustering properties of galaxies (e.g., 
Berhnd & Weinberg 2002; Berlind et al. 2003). 

No such non-linear relation is evident between 
A^200 and -^200 7 the Ks band luminosity inside 
i?200- This result shows that the non-linearity 
of the halo occupation function is not driven by 
variations in the luminosity function. 

More massive virialized halos have larger mass- 
to-light ratios. This result follows logically from 
the two prior points. Our M/L — M relation 
agrees with previous determinations (Bahcall 
& Comerford 2002, L03). These results signify 
that the efficiency of galaxy formation decreases 
(and/or that the efficiency of galaxy disruption 
increases) with increasing halo mass and/or virial 
temperature. 



• We investigate possible systematic effects and 
conclude that dark matter is more concentrated 
than stellar mass contained in galaxies. This 
result could arise either from different clustering 
properties of dark matter and baryonic matter 
or from variations in the efficiency of converting 
baryonic matter into galaxies. The cluster environ- 
ment seems to be either less efficient at converting 
baryons into galaxies or more efficient at disrupting 
galaxies than less dense environments. Such a 
difference is predicted by simulations of ACDM 
cosmologies where processes such as tidal stripping 
and dynamical friction disrupt galaxies in clusters 
(Kravtsov & Klypin 1999; Colin et al. 1999) and 
supported by observations of significant numbers 
of intergalactic stars in clusters. Alternatively, 
the heating of the intracluster medium may cut 
off the supply of cold material needed to form 
stars (e.g., Blanton et al. 1999; Balogh et al. 2000, 
and references therein), thus lowering the star 
formation efficiency in cluster galaxies. 

• Assuming the mass-to-light ratios at large radii 

are similar to the global value, we estimate 
^Im = 0.10±0.03 (1-<J statistical uncertainty) using 
the SDSS luminosity density with appropriate 
color corrections or flm = 0.13 ± 0.03 (l-cr 
statistical uncertainty) from the 2dFGRS. We 
suggest that the 2dFGRS and the CfA/SSRS2 
surveys sample local underdensities. Uncertainties 
in the luminosity density, especially at infrared 
wavelengths, contribute a significant amount of the 
systematic uncertainty in estimating Q^n- These 
estiniatcis of n„, are small compared with other 
recent estimates from the microwave background, 
the galaxy power spectrum, and supernovae. 
However, they agree well with other estimates 
based on cluster mass-to-light ratios (Carlberg 
et al. 1996, 1997; BahcaU et al. 2000; Girardi 
et al. 2000; Bahcall & Comerford 2002), cluster 
abundances (Reiprich & Bohringer 2002; Bahcall 
et al. 2003, but see Scliuecker et al. 2003) and weak 
lensing (Kaiser et al. 2004; Wilson et al. 2001; 
Hoekstra et al. 2001; Gray et al. 2002). We discuss 
possible systematic effects that could cause our 
result to be anomalously low. Reconciling these 
estimates of il„j by invoking bias requires that the 
typical value of M/ Lk at the smallest densities 
we probe w 3/9c is a factor of 2-3 smaller than the 
global value. For instance, if galaxy formation 
occurs nearly exclusively above a density threshold 
5 ~ 10, the mass-to-light ratios in cluster outskirts 
may underestimate the global value. 

One promising future direction is to study clusters at 
moderate redshifts where weak lensing provides an inde- 
pendent mass estimate (Kneib et al. 2003). Comparing 
lensing mass profiles to caustic mass profiles will constrain 
unknown systematics in both techniques. For instance, a 
sheet of mass of uniform density produces no lensing signal 
(the mass-sheet degeneracy), but this mass should be evi- 
dent in the galaxy kinematics. Conversely, foreground and 
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background structures may produce a weak lensing signal 
but would not affect the kinematics of the infall region. 

Infall regions are interesting environments for studying 
the evolution of galaxy populations. We show here that in- 
fall regions are important in constraining models of galaxy 
bias or antibias and fi^. If other methods yield a precise 
measurement of flm, the changes in mass-to-light ratios 
with environment provide important clues to the forma- 
tion and evolution of galaxies. 
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Table 4 

CAIRNS Near-Infrared Spectroscopic Completeness 



Cluster 


r200 


rt 






fnoz 




KEkM 








h-^Mpc 


h-lMpc 


/i-iMpc 


mag 




mag 


mag 


mag 




A119 


1.07 


5.4 


5.4 


12.9 


0.096 


0.014 


-0.125 


-22.70 


0.616 


A168 


1.09 


5.5 


5.5 


13.1 


0.092 


0.013 


-0.127 


-22.50 


0.664 


A496 


0.98 


4.2 


4.0 


12.6 


0.042 


0.050 


-0.094 


-22.38 


0.690 


A539 


1.03 


4.3 


2.5 


12.4 


0.086 


0.062 


-0.083 


-22.31 


0.704 


A576 


1.42 


6.0 


4.3 


12.9 


0.036 


0.028 


-0.110 


-22.38 


0.690 


A1367 


1.18 


5.2 


5.2 


12.3 


0.033 


0.008 


-0.062 


-21.75 


0.801 


Coma 


1.50 


7.4 


7.4 


12.7 


0.028 


0.003 


-0.068 


-21.61 


0.820 


A2199 


1.12 


5.3 


4.1 


13.0 


0.027 


0.004 


-0.087 


-21.78 


0.800 


A194 


0.69 


3.3 


3.3 


12.2 


0.044 


0.015 


-0.052 


-21.42 


0.844 



Assuming a luminosity function with = —23.77 and a = —1.14 as in the 2dF/2MASS luminosity function for extrapolated Kg magnitudes. 
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Fig. 1. — Completeness of the CAIRNS spectroscopic catalogs versus absolute magnitude Mk^ (falling thick solid lines and scales on left). 
Vertical bars indicate the spectroscopic completeness limits. The scales on the right show the number of galaxies per 0.2 magnitude bin (rising 
thin solid lines). The dashed lines show the number of galaxies with redsliifts and the dash-dot lines show the number of cluster galaxies with 
the upper and lower lines indicating the maximum and minimum number of members. The scales are identical in all panels. Clusters are 
ordered left to right, top to bottom, in decreasing X-ray temperature. 
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Fig. 2. — Galaxy luminosity functions for all galaxies within the infall regions of the CAIRNS clusters in Kg band. The solid line is the 
field LF for comparison with arbitrary absolute normalization but relative normalization scaled by Ntot, the number of galaxies brighter than 
= —22.77 within the limiting radius of the caustics. 
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Fig. 3. — Same as Figure 2 for galaxies within il200- Relative normalizations are made using N200, the number of bright galaxies projected 
within il200- 
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Fig. 4. — Same as Figure 2 for galaxies outside il200- Relative normalizations are made using ATj^y, the number of bright galaxies projected 
outside R200- 




Fk;. 5. Combined ifs-band luminosity functions for the CAIRNS clusters. Open circles show the total LF, squares show the LF inside 
^200 1 and triangles show the LF outside i?2()0- Solid lines show the shape of the 2dF/2MASS A's-band LF with arbitrary normalization. The 
LFs inside and outside iJ200 offset for clarity. The vertical line indicates + 1, the completeness limit of the survey. 
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Fig. 6. — Ratio of the combined CAIRNS JsTg-band luminosity function of galaxies projected outside R200 to that of galaxies projected 
inside ii200 ■ Upper limits show the ratio if one galaxy were present outside il200 • The vertical solid line shows the approximate limit of the 
survey Mjc^ = —22.77, one magnitude fainter than for field galaxies. The horizontal line shows the ratio of all galaxies (M^, < —22.77) 
projected outside R200 to all galaxies {Mk, < —22.77) inside R200- Errorbars indicate l-cr Poissonian uncertainties. 
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Table 5 

Luminosity Function Parameters 



LF 


Filter 


M* 


a 


CAIRNS Virial 


Ks 




1 Qcr+U.ib 
i.Od_g jg 


CAIRNS Infall 


Ks 


-23.86to2? 


1 90+0.17 

-'-■''"^-0.15 


CAIRNS Total 


Ks 


-23.97l°;ii 


-1 26+°-°^ 


2dF/2MASS Extrap 
2dF/2MASS Kroni 
CfA/2MASS2 


Ks 
Ks 
Ks 


-23.77 ± 0.03 
-23.64 ± 0.03 
-23.63 ± 0.05 


-1.14 ±0.05 
-0.96 ± 0.05 
-1.09 ± 0.06 


CAIRNS Virial 


J 


„,.+U.iy 

zo.zo_Q 23 


-1 26+"-^^ 


CAIRNS Infall 


J 


—99 70+01''' 


, 1 9+0.16 


CAIRNS Total 
2dF/2MASS Extrap 
2dF/2MASS Kron^ 


J 
J 
J 


-23.00t»:l^ 
-22.70 ± 0.02 
-22.56 ±0.02 


-1.07 ±0.03 
-0.93 ± 0.04 



^Wc shift the Mj^ of COl by -0.20 to convert from Kron magnitudes (from J band Kron magnitudes and J — Kg colors) to total (extrapolated) 

2MASS Ks magnitudes (see Figure 5 in COl). 
^Wc shift the AfJ^ of KOI by -0.05 to convert from isophotal to Kron magnitudes and then by -0.20 to convert from COl Kron magnitudes to 

2MASS Ks total magnitudes. 
^We shift the Mj of COl by -0.20 to convert from Kron magnitudes to total (extrapolated) 2MASS J magnitudes. 



Table 6 

CAIRNS Rkilxess 



Cluster 




'f'rii a X 

ft-iMpc 


-L200 ^ 




iV200'' 




A119 


1.07 


5.4 


6.66 


13.75 


47 


69 


A168 


1.09 


5.5 


3.61 


11.69 


35 


83 


A496 


0.98 


4.0 


3.94 


7.86 


38 


40 


A539 


1.03 


2.5 


3.32 


5.43 


31 


21 


A576 


1.42 


4.3 


6.39 


12.35 


58 


59 


A1367 


1.18 


5.2 


4.44 


9.75 


44 


51 


Coma 


1.50 


7.4 


10.40 


21.96 


89 


107 


A2199 


1.12 


4.1 


5.15 


17.43 


52 


112 


A194 


0.69 


3.3 


1.72 


4.02 


16 


20 



^Luminosities are corrected for faint galaxies by dividing by from Table 4. 

^Richness is for galaxies brighter than Mjc^ = —22.77 + 51og?i, equivalent to one magnitude fainter than for field galaxies. 



Table 7 
CAIRNS Scale Radii 



Cluster 


^200 


^max 


0.M 










a.L 








NFW 


NFW 


NFW 


Hern 


Hern 


Hern 


A119 


1.07 


5.4 


0.17 


0.33 


0.19 


0.58 


1.24 


1.00 


A168 


1.09 


5.5 


0.21 


0.63 


0.75 


0.65 


2.17 


2.50 


A496 


0.98 


4.0 


0.07 


0.22 


0.13 


0.31 


0.81 


0.75 


A539 


1.03 


2.5 


0.07 


0.08 


0.07 


0.25 


0.49 


0.45 


A576 


1.42 


4.3 


0.13 


0.35 


0.25 


0.43 


1.22 


1.12 


A1367 


1.18 


5.2 


0.07 


0.27 


0.33 


0.29 


1.18 


1.30 


Coma 


1.50 


7.1 


0.15 


0.35 


0.28 


0.50 


1.10 


i.:^2 


A2199 


1.12 


i.l 


0.1 .J 


0.7.5 


0.90 


0. 17 


2.13 


2.58 


A194 


0.69 


3.3 


0.11 


0.20 


0.27 


0.35 


0.93 


1.14 
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Table 8 

CAIRNS Mass-to-Light Ratios 



Cluster 


'"200 

fe-iMpc 


h-lMpc 


^-200 
10l2h-2LQ 


Ltot 

lOl2h-2L0 


(M/Lnjtot 
HMq/Lq 


(M/Lx J200 

hMQ/LQ 


(M/Litjinf 

HMq/Lq 


A119 


1.07 


5.4 


6.66 


13.75 


46±27 


43±18 


49±21 


A168 


1.09 


5.5 


3.61 


11.69 


56±16 


83±21 


44±10 


A496 


0.98 


4.0 


3.94 


7.86 


38±14 


56±15 


21±6 


A539 


1.03 


2.5 


3.32 5.43 


59±25 


76±31 


31±9 




A576 


1.42 


4.3 


6.39 


12.35 


71±14 


104±19 


36±7 


A1367 


1.18 


5.2 


4.44 


9.75 


57±18 


86±24 


34±8 


Coma 


1.50 


7.4 


10.40 


21.96 


75±21 


75±12 


75±20 


A2199 


1.12 


4.1 


5.15 


17.43 


33±11 


63±17 


21±5 


A194 


0.69 


3.3 


1.72 


4.02 


37±28 


44±24 


32±19 



Table 9 
Estimates of 



Technique 




Reference 


CAIRNS Total M/L + 2dF 


0.18 ±0.04 




CAIRNS Virial M/L + 2dF 


0.24 ±0.05 




CAIRNS Infall M/L + 2dF 


0.13 ±0.03 




CAIRNS Total M/L + SDSS 


0.14 ±0.05 




CAIRNS Virial M/L + SDSS 


0.18 ±0.06 




CAIRNS Infall M/L + SDSS 


0.10 ±0.03 




WMAP + 2dF 


0.27 ±0.04 


1 


WMAP + SDSS 


0.30 ± 0.04 


2 


WMAP + Other CMB 


0.1 - 0.5 (95%) 


3 


Type la SNe 


0.28 ±0.05 


4 


Type la SNe 


0.25+Q og(stat) ±0.04(sys) 


5 


Cluster Abundance 


12+°'"' 


6 


Cluster Abundance 


0.17 ±0.05 


7,8 


Clus. Abun. +Clustering 


0.34 ± 0.03(stat) ± 0.09(sys) 


9 


Weak Lcnsing (Groups) 


0.13 ±0.07 


10 


Weak Lcnsing (Superclusters) 


<o.i 


11,12 


Weak Lensing (Blank) 


K» 0.1 


13 


Gas Mass Fraction 


0.28 ±0.03 


14,15 


CMB + Power Spectrum + BBN + fg 


0.33 ± 0.04 


16 



References. — (1) Spergel et al. (2003); (2) Tegmark et al. (2003); (3) Bridle et al. (2003); (4) Tonry et al. (2003); (5) Knop et al. (2003); 
(6) Reiprich & Bohringer (2002); (7) Bahcall et al. (2003); (8) Bahcall & Bode (2003); (9) Schuecker et al. (2003); (10) Hoekstra et al. (2001); 
(11) Kaiser et al. (2004); (12) Gray et al. (2002); (13) Wilson et al. (2001); (14) Ettori & Fabian (1999); (15) L03; (16)Turner (2002). 



28 



-25 -24 -23 -22 -21 -25 -24 -23 -22 -21 -25 -24 -23 -22 -21 




-25 -24 -23 -22 -21 -25 -24 -23 -22 -21 -25 -24 -23 -22 -21 

Mj -H 5 log h 



Fig. 7. — Galaxy luminosity functions for all galaxies within the infall regions of the CAIRNS clusters in J band. The solid line is the 
field LF for comparison with arbitrary absolute iiorinalization but relative normalization scaled by Ntot, the number of galaxies brighter than 
Mj = —21.70 within the limiting radius of the caustics. 
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Fig. 8. — Same as Figure 7 for galaxies within i?200- 
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Fig. 9. — Same as Figure 7 for galaxies outside R200- 
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Fig. 10. — Combined J-band luminosity functions for the CAIRNS elusters. Open circles show the total LF, squares show the LF inside 
i?200i a-nd triangles show the LF outside -R200. Solid lines show the shape of the 2dF/2MASS J-band LF with arbitrary normalization. The 
LFs inside and outside i?200 offset for clarity. The vertical line indicates Mj + 1, the approximate completeness limit of the survey. 
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Fig. 11. — Ratio of the combined CAIRNS J-band luminosity function of galaxies projected outside Rioo to that of galaxies projected 
inside il200 . Upper limits show the ratio if one galaxy were present outside il200 . The vertical solid line shows the approximate limit of the 
survey Mj = —21.7, one magnitude fainter than Mj for field galaxies. The horizontal line shows the ratio of all galaxies {Mj < —21.7) 
projected outside R200 to all galaxies {Mj < —21.7) inside R200- Errorbars indicate l-cr Poissonian uncertainties. 
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Fig. 12. — Absolute magnitude versus (projected) clustercentric distance in units of r2oo- Dotted lines indicate for field galcixies and 

solid lines indicate the spectroscopic completeness limits. 
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Fig. 13. — A576 R — Ks color versus projected distance from the cluster. The symbols show galaxies of different luminosities. Large 
filled squares, small filled squares, hexagons, and stars represent galaxies in the magnitude bins Mpc^ < —23.77, —23.77 < < —22.77, 

-22.77 < Mx, < -21.77, and -21.77 < Mk. < -20.77 respectively. Note that the latter two bins are not complete. 
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(mag) 

Fig. 14. — A576 R — Ks color (in a circular aperture of radius 15") versus absolute Ks band magnitude within this aperture. Squares are 
galaxies inside i?200 a-nd crosses arc galaxies outside i?200 ■ The vertical line shows the magnitude limit of the 2MASS spectroscopic catalog 
while the slanted line shows the limit of the R band spectroscopic catalog. 
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Fig. 15. — Near-infrared J — Ks color versus projected distance from the cluster. Circles and crosses represent galaxies in the magnitude 
bins Mk, < -23.77 and -23.77 < Mk, < -22.77 respectively. 
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Fig. 16. — Normalized histograms of near-infrared J — Kg colors inside (solid lines) and outside (dash-dotted lines) ii200. 
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Fig. 17.- The ncai-infrarcd color-magnitude relation for the CAIRNS clusters. Crosses are galaxies projected inside R200 a^nd triangles 
arc those projected outside i?20()- The solid and dotted vertical lines indicate the spectroscopic completeness limits and . We overplot a 

fiducial color-magnitude relation with slope — O.Olmag mag~^. 
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Fig. 18. — Surface number density profile of cluster members brighter than Mn^ = —22.77 + Slogft, equivalent to + 1 for field galaxies. 
The dashed line shows a singular isothermal sphere. The solid and dash-dotted lines show the best-fit NFW and Hernquist profiles. 
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Fig. 19. — Surface luminosity density profiles for cluster members brighter than Mk^ = —22.77 + 51og/i, equivalent to + 1 for field 

galaxies. The dashed line shows a singular isothermal sphere. The solid and dash-dotted lines show the best-fit NFW and Hernquist profiles. 
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Fig. 20. — Mass-to-light ratio in Ks band as a function of radius for the CAIRNS clusters. The solid lines show the caustic mass profile 
M(< r) divided by the projected luminosity profile (< Rp) ^^nd the shaded regions indicate the associated l-cr uncertainties. The open 
squares show the mass-to-light ratio M{r,r + dr)/ Lji^{Rp,Rp + dRp) in radial shells. The dash-dotted line shows the projected best-fit 
Hernquist mass profile Mh{< Rp) divided by (< Rp)- The stars show the mass-to-light ratio calculated using the X-ray temperature 
and the mass-temperature relation to estimate the mass. The lower point shows Mx{< J'50o)/i'ifa (< R500) and the upper point shows 
Mx{< R50o)/Lks{< R500) assuming Mx{< R500) = l.SMjf (< rsoo) as is true for an NFW mass profile with c=5. 
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Fig. 21. — Number of bright galaxies N2oo,P projected within il200 versus M200 for the CAIRNS clusters. The solid line shows the bisector 
of the two ordinary least-squares fits to the data. The dashed line shows N200 oc M200 . 
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Fig. 22. — Mass to light ratio versus M200 for the CAIRNS clusters. The solid curve shows the relation for X-ray clusters in LOS (converted 
to M20o)i the dotted line the relation for 2MASS clusters with the model of Kochanck et al. (2003), and the dash-dotted line shows the 
M/L — Tx relation of Bahcall & Comerford (2002) converted to M/L — M using the mass-temperature relation of Finoguenov et al. (2001). 
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Fig. 23. — Luminosity within il200 versus number of bright galaxies within ii200. The soUd Une shows the bisector of the two ordinary 
least-squares fits and the dashed Une shows iV2oo <x. L200 ■ 
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Fig. 24. — Velocity dispersion profiles for the CAIRNS clusters tciken from Paper I. The bars on the abscissa indicate r2oo and ft- The 
dash-dotted lines are the VDPs of the Hernquist mass profiles that best fit the caustic mass profiles assuming isotropic orbits. The solid lines 
show the Hernquist VDP predicted from the surface number density profiles assuming isotropic orbits and a constant ratio of mass to number 
density. 
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Fig. 25. — Comparison of Ks band and R band mass-to-light profiles for A576. The vertical scale on the right is for R band. The Kg 
band mass-to-light profile is the lower set of solid lines with l-cr uncertainties. The dash-dotted line and shaded region show the R band 
mass-to-light profile and l-cr uncertainties from Rines et al. (2000). The dashed line shows the R band mass-to-light profile and surrounding 
solid lines show the l-cr uncertainties calculated from bright galaxies assuming a constant fraction of light in fainter galaxies. 



